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Koordinatensysteme in Feuerleitgeraten 


Von Amsros P, SpetseEr, Ziirich?) 


1. Einleitung 


Die vorliegende Arbeit behandelt Feuerleitgerate fiir die terrestrische Flug- 
zeugabwehr. Das Hauptproblem bei diesen Geraten ist bekanntlich die Extra- 


polation, das heisst die Aufgabe, den Standort des Ziels auf langere Zeiten (zum ~~ 


Beispiel 20s) vorauszuberechnen. Im Folgenden wird die Wahl geeigneter - 
Koordinatensysteme fiir die verschiedenen Teile eines Feuerleitgerates be- - 
handelt. Es sind nur die wichtigeren Schritte des ganzen Rechenprozesses 
beriicksichtigt. Ferner wird nur die lineare, nicht aber die kurvenférmige ~ 
Extrapolation behandelt. 


2. Verwendete Bezeichnungen 


(Vektoren sind durch fette Buchstaben gekennzeichnet.) 


Messpunkt 


M 
= 
y |  Kartesische Messpunktkoordinaten 
ae 
a Seite 

e 


Entfernung 


ey, Karten-(Horizontal-) Entfernung fee he 
M Ortsvektor 
v Geschwindigkeitsvektor 
v;, Horizontalgeschwindigkeit 
v* Geglattete Geschwindigkeit 
x Kurswinkel 
pb Flugrichtung 
: se Treffpunkt 
; Ty | : ; 
: Yr | Kartesische Treffpunktkoordinaten 
hr 
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xr Seite 1 
Ar Lagewinkel | 
er Entfernung des Treffpunktes 
Cr Karten-(Horizontal-) Entfernung | 
Ay Ortsvektor 

E Elevation 

o Schusswinkel 

a Flugzeit 

oO Tempierung 

Uo Geschossanfangsgeschwindigkeit 

C) Luftgewicht 

Rk Ladeverzug 

P Parallaxvektor 

iL Linearer Differentialoperator 


Figur 1 Figur 2 
Erlauterung der Bezeichnungen: Horizontal- Erlauterung der Bezeichnungen: Visierstrahl 
projektion des Vorhaltedreieckes. und Flugbahn in der Schussebene. 


Die Figuren 1 und 2 dienen zur Erlauterung einzelner dieser Bezeichnungen. 


3. Die Elemente eines Feuerleitgerdtes 


Bei der Wahl von geeigneten Koordinatensystemen kann der Konstrukteur 
nicht jede mathematische Operation fiir sich getrennt behandeln, sondern er 
muss das Feuerleitgerat als ganzes betrachten und muss darnach streben, 
zwischen Materialaufwand einerseits und erreichter Genauigkeit andererseits 
ein méglichst giinstiges Verhaltnis zu erzielen. Es ist von folgenden Randbe- 
dingungen auszugehen: Das Visier (gleichgiiltig, ob es sich um ein optisches 
Visier oder um Radar handelt) soll in Azimutalkoordinaten «, A, e arbeiten; 
dem Geschiitz miissen (abgesehen von Ausnahmefallen) die Koordinaten ar, 
€ sowie die Tempierung # zugefiihrt werden. Der Grund, warum zwischen 
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diesen beiden Koordinatengruppen uberhaupt eine Transformation nétig ist 
hegt darin, dass der Messpunktvektor und der Vorhaltevektor addiert werden 
miissen. Bekanntlich ist die Addition zweier Vektoren nur in kartesischen 
Koordinaten komponentenweise méglich. Es muss also entweder eine Trans- 
formation in kartesische Koordinaten vorgenommen werden, oder es miissen 
die wesentlich komplizierteren Additionsformeln fiir Vektoren in krummlinigen 
Systemen zur Verwendung gelangen [1]?) 


2, Koord. 


M=y — y=L(v*) v™r 


Figur 3 


Blockschema der Elemente eines Feuerleitgerates (ohne Korrekturteil). 


M/2,i,e) 


Figur 4 


Variante von Figur 3: Die Koordinatentransformatoren befinden sich nur im Geschwindigkeits- 
zweig. 


Figur 3 veranschaulicht den prinzipiellen Aufbau eines Feuerleitgerates. 
Fiir jede Funktion ist ein gesonderter Block gezeichnet. Die in der Praxis 
verwendeten technischen Lésungen gestatten oft die Zusammenfassung von 
zwei oder mehreren der aufgefiihrten Funktionen in einem einzigen Apparat. 
Fiir die rechnerische Behandlung ist aber festzuhalten, dass jeder Block einen 
gesonderten mathematischen Prozess darstellt. 

Das Visier liefert den Messpunktvektor M. Dieser wird einer Koordinaten- 
transformation unterzogen und anschliessend differenziert, indem das zeit- 
liche Differential gleich dem Geschwindigkeitsvektor v ist. Dieser Geschwin- 
digkeitsvektor wird geglaittet, damit kleine Stérungen sich nicht bemerkbar 
machen, und im Multiplikator mit der Geschossflugzeit t+ miultipliziert. Als 
Resultat entsteht der Vorhaltevektor. Nach folgender Addition entsteht als- 

dann der Treffpunktvektor T: 


T= M+ 0* x. 


2) Die Ziffern in eckigen Klammern verweisen aut das Literaturverzeichnis, Seite 16. 
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Eine neue Koordinatentransformation fiihrt den entstandenen Vektor in eine 
Form iiber, die zur Verarbeitung im ballistischen Auswerter geeignet ist. Der 
ballistische Auswerter liefert die Geschiitzkoordinaten. 


4. Glattung der Geschwindigkeit 


Bevor wir die Koordinatensysteme betrachten, in denen der Geschwindig- 
keitsvektor dargestellt wird, muss das Prinzip der Glattung erlautert werden, 
da Bildung und Glattung von v eng miteinander zusammenhangen. Die Glat- 
tung verfolgt den Zweck, so weit als méglich das im Messpunktvektor enthal- 
tene Stérsignal vom Nutzsignal zu trennen. Das Stérsignal rihrt von unge- 
nauer Arbeit des Visiers (eventuell des Richtpersonals) und der nachfolgenden 
Rechenelemente her. Nun ist das Problem der Trennung von Nutz- und St6r- 
signal ausfiihrlich behandelt worden [2], unter der Voraussetzung, dass man 
sich auf die Verwendung von linearen Differentialoperatoren, dargestellt zum 
Beispiel durch lineare elektrische Filter, beschrankt. Nach dieser Theorie 
miissen sowohl vom Nutz- als auch vom Stoérsignal diejenigen statistischen 
Eigenschaften bekannt sein, welche in den Korrelationsfunktionen enthalten 
sind. Daraus wird ersichtlich, dass die Arbeit eines Feuerleitgerates als stati- 
scher Prozess aufzufassen ist, worauf weiter unten eingegangen werden soll. 
Die statistischen Eigenschaften des St6rsignals, insoweit sie durch die Korre- 
lationsfunktion gekennzeichnet sind, kénnen als zeitlich konstant betrachtet 
werden, da die Entstehung der Stérungen immer auf die gleiche Ursache 
_ zuriickzufiihren ist. Dagegen ist die Korrelationsfunktion des Nutzsignales, das 
heisst des Messpunktvektors, weiten Schwankungen unterworfen. Ihre Be- 
schaffenheit hangt ab von der Art des zu bekampfenden Zieles, von dessen 
Geschwindigkeit, von der momentanen taktischen Lage, ferner vom Verhalten 
des Feuerleitenden, indem es nicht auf dasselbe herauskommt, ob ein Ziel im 
Verlaufe eines Vorbeifluges wahrend kiirzerer oder wahrend langerer Zeit 
bekampft wird. Da nun die Eigenschaften des optimalen Filters von den Korre- 
lationsfunktionen abhangen, ist es auf keinen Fall angangig, in einem Feuer- 
leitgerat diese Eigenschaften in Form von unveranderlichen Konstanten ein- 
zubauen. 

Das heute vorwiegend verwendete Filter bildet die geglattete Geschwindig- 
jkeit D* nach folgender Formel: 
v= Liv*). (1) 


L kennzeichnet hier einen linearen Differentialoperator. Am haufigsten setzt 
man: a 


Dieser Operator bewirkt, dass v* sich zusammensetzt aus dem momentanen 
Wert von U sowie aus gewogenen Anteilen der Geschwindigkeiten aller ver- 
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gangenen Zeiten, wobei die Gewichte um so kleiner werden, je weiter der be- 
treffende Wert zuriickliegt. Die Zeitkonstante 7y-ist ein Mass dafiir, wie weit 
in die Vergangenheit zuriickgegriffen wird. Es bedeutet einen wesentlichen 
Vorteil fiir ein Feuerleitgerat, wenn die Grésse von JT, durch den Feuerleitenden 
eingestellt werden kann. Die Ansichten iiber die beste Grésse von Ty fiir ver- 
schiedene Verhaltnisse gehen sehr stark auseinander, und es werden Werte 
von 0,5 s bis 30 s vertreten. Die Beurteilung dieser Auffassungen erfordert eine 
eingehende Untersuchung der verschiedenen damit zusammenhingenden Vor- 
gange, doch kann allgemein gesagt werden, dass bei schnell bewegten Zielen 
und gut arbeitendem Visier kleine Zeitkonstanten zu wahlen sind, wahrend bei 
langsam bewegten Zielen, die eventuell noch Ausweichmanéver durchfiihren, 
bei unregelmassig arbeitendem Visier und schliesslich bei kurzzeitigem Ver- 
schwinden des Zieles, etwa hinter Wolken, gréssere Zeitkonstanten die besten 
Resultate ergeben. 

Die in [2] dargelegte Theorie gestattet selbstverstandlich auch die Behand- 
lung komplizierterer Operatoren als des durch (2) beschriebenen, jedoch ist 
sie auf lineare Operatoren beschrankt. Es ist aber nicht daran zu zweifeln, 
dass beste Resultate an Feuerleitgeraten nur mit nichtlinearen Filtern erreicht 
werden kénnen. Daher ist es unwahrscheinlich, dass mit den heute verfiig- 
baren mathematischen Methoden auf wissenschaftlicher Grundlage ein wirk- 
lich optimales Filter gefunden werden kann. Man ist vielmehr auf Lésungen 
angewiesen, deren Wirksamkeit empirisch tiberpriift wird. Eine Annaherung 
an ein Filter von sehr allgemeiner Natur wird schon seit langer Zeit in vielen 
Feuerleitgeraten fiir gewisse Komponenten von © verwendet, namlich das 
Nachstellen eines Folgezeigers durch einen Bedienungsmann. Dieser Bedie- 
nungsmann kann fiir seine spezielle Aufgabe geschult werden und kann ausser- 
dem von Fall zu Fall durch den Feuerleitenden itber Art des zu bekampfenden 
Zieles sowie die taktische Lage und deren mutmassliche Entwicklung orientiert 
werden und entspricht daher einem nichtlinearen Filter mit zeitlich verander- 
lichen Eigenschaften, dessen Vollkommenheit durch eine maschinelle Vorrich- 
tung wohl nur schwer zu erreichen ist. 

Die Verwirklichung eines manuellen oder automatischen Filters stésst auf 
gewisse praktische Schwierigkeiten, da die anzuwendenden Operatoren in der 
Regel nur fiir skalare Gréssen, nicht aber fiir Vektoren dargestellt werden 
k6nnen. Man muss daher auch fiir v ein Koordinatensystem wahlen. Hier ist 
es nun ein grosser Vorteil, wenn ein System gefunden werden kann, in welchem 
die Komponenten von v bei geradlinigem Flug mit konstanter Geschwindigkeit 
konstant sind. Dies kann wie folgt begriindet werden. Durch Anwendung des 
Differentialoperators L soll ein Signal, welches von keinerlei Storungen behaftet 
ist, nicht verandert werden. Dies rithrt von der elementaren Forderung her, die 
an ein Filter dieser Art gestellt werden kann, dass gute Signale nicht verschlech- 
tert werden diirfen. Nun lassen lineare Differentialoperatoren, insoweit sie 
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hier tiberhaupt in Betracht kommen, ein Signal nur dann unverandert, wenn 
es gleich einer Konstanten ist. Daraus resultiert die Bedingung fiir v, dass die 
Darstellung in einem raumfesten Koordinatensystem erfolgen muss. Bewegte 
Systeme ergeben variable Komponenten, und die Anwendung eines Filters auf 
solche Komponenten wiirde einer Verschlechterung von guten Signalen gleich- 
kommen. Dementsprechend besitzen Gerate, welche v in ortsfesten Koordi- 
naten darstellen, einen wesentlichen Vorsprung. 

Die statistische Theorie in [2] beschrankt sich tibrigens in ihrer Anwendung 
nicht auf die Glattung von v. Vielmehr ist es méglich, als Aufgabe fiir das 
Filter die direkte Bildung von T aus M zu formulieren, so dass ein einziger 


Apparat gleichzeitig eine Extrapolation und eine Glattung ausfihrt. Wegen — 


der erwaihnten Beschrankung auf lineare und zeitlich unabhangige Vorrich- 


tungen ist es aber unwahrscheinlich, dass das Problem in dieser allgemeinen ~ 


Form eine brauchbare Lésung findet. 


5. Bildung der Geschwindigkeit 


a) Transformation des Messpunktvektors in kartesische Koordinaten 


Die mathematisch einfachste Lésung fiir die Bildung von v ist eine Trans- 
formation des Messpunktvektors M in kartesische Koordinaten x, y, 4. Durch 
Differentiation jeder einzelnen Komponente entsteht sofort v. 

Diese Differentiation wird haufig mit Reibradgetrieben durchgefiihrt. 
Reibradgetriebe sind an sich integrierende Rechengeradte, werden aber durch 
eine Art Gegenkopplungsschaltung als Differentiatoren verwendet. Sie bilden 
dann nicht das exakte Differential, sondern eine Funktion, die sich vom Dif- 
ferential gerade durch den Operator (2) unterscheidet. Sie stellen also Differen- 
tiatoren und Filter zugleich dar. Im allgemeinen lasst sich aber Ty hier nicht 
verandern, so dass der Wert dieser Filterung stark herabgemindert ist. Andere 
Gerate bilden den Differentialquotienten in der Weise, dass sie die mittlere Ge- 
schwindigkeit iiber ein Zeitintervall von zum Beispiel 2s verwenden. Dies 
geschieht mit mechanischen Anordnungen, welche als Stop-Tachometer be- 
zeichnet werden. Die entstehende Funktion entspricht der Anwendung eines 
zeitabhangigen Differentialoperators, der eine gewisse Ahnlichkeit mit (2) 
besitzt, wobei T, ungefahr dem Messintervall entspricht. Die so entstehenden 
Geschwindigkeiten haben aber den Nachteil, dass sie sich sprungweise andern. 


b) Geometrische Abbildung des Geschwindigkeitsvektors 


Nachfolgend wird die Darstellung von v in raumfesten Zylinderkoordinaten 
mit den Komponenten v,, x, v, = h beschrieben (fiir die Erklarung des Kurs- | 


winkels % siehe Figur 1). Die drei Komponenten werden aut ganz verschiedene 
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Art gebildet. Zur Ermittluing der Horizontalgeschwindigkeit v, dient eine 
Schlepprolle gemiss Figur 5. Diese wird nach Massgabe der Horizontalprojek- 
tion des Flugweges tiber eine raumfeste Platte gefithrt. Durch die Eigenschaften 
der Schlepprolle stellt sich an ihrer vertikalen Achse der Kurswinkel x ein, und 
zwar wiederum nicht exakt, sondern unter Beriicksichtigung eines Operators 
ahnlich demjenigen in (2), wobei jedoch 7, vom Betrag der Horizontalgeschwin- 
digkeit abhangt. Dieser Betrag ist durch die Umdrehungsgeschwindigkeit der 
Rolle dargestellt und wird dadurch gemessen, dass das Gerat die Zeit zwischen 
der Schliessung zweier aufeinanderfolgender Kontakte, die entlang dem Umfang 


Figur 5 


Schlepprolle, drehbar um eine horizontale und eine vertikale Achse. Die fiinf angedeuteten Kon- 
takte dienen zur Messung der Geschwindigkeit. 


der Rolle angebracht sind, auswertet. Das Verfahren ist ahnlich einem Stop- 
Tachometer, jedoch mit konstanter Wegstrecke, nicht mit konstanten Zeit- 
intervallen. Das ermittelte v, wird sich also sprunghaft andern, und hier kann 
nun das Nachstellen eines Folgezeigers durch einen Bedienungsmann einge- 
schaltet werden, was den im letzten Abschnitt geschilderten grossen Vorteil 
bringt. 

Es sei noch erwahnt, dass man gelegentlich als Unterlage fiir die Schlepp- 
rolle an Stelle einer Platte eine Kugel wahlt. Die vertikale Achse der Schlepp- 
rolle ist dann ortsfest, und die Kugel wird mit Hilfe angepresster Reibrader 
nach Massgabe des Flugweges rotiert. Diese Anordnung beansprucht bedeutend 
weniger Platz. ; 

Die noch fehlende Komponente v, = / kann auf irgendeine der im Abschnitt 
5a beschriebenen Arten ermittelt werden. 


c) Bildung von v in rechtwinkligen bewegten Koordinaten © 


Als Koordinatensystem fiir v kann ein rechtwinkliges bewegtes System 
-gewahlt werden, dessen eine Achse mit der Visierlinie oder mit deren Horizon- 
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talprojektion zusammenfallt, wahrend die zweite Achse horizontal liegt. Dem- 
entsprechend steht die dritte Achse senkrecht zur Visierlinie bzw. vertikal. 
Die Bildung der Komponenten von Uv erfolgt aus den Koordinaten «, A, e des 
Messpunktes sowie aus deren zeitlichen Ableitungen. Dieser Vorgang ist ein- 
fach durchzufiihren, besonders deshalb, da die zeitlichen Differentiale der Visier- 
werte ohnehin verfiigbar sind, falls mit einer Weg-Geschwindigkeits-Steuerung 
gearbeitet wird. Eine solche Steuerung liefert die zeitlichen Ableitungen mit 
guter Genauigkeit und mit einem Minimum an Stérungen. Diesem bedeutenden 
Vorteil steht der Nachteil gegeniiber, dass die Komponenten der Geschwindig- 
keit auch bei geradlinigem Flug nicht konstant sind. Dies verschlechtert die 
Wirkung der Glattung, wie im Abschnitt 4 dargelegt wurde. 


d) Transformation des Geschwindigkeitsvektors in ortsfeste kartesische Koordinaten 


Ein weiteres Verfahren der Darstellung von v besteht darin, diesen Vektor 
zunachst in einem bewegten rechtwinkligen Koordinatensystem durch Bildung 
der Komponenten é, e« cos, edz erzeugen. Alsdann wird der entstandene 
Vektor in raumfeste rechtwinklige Koordinaten x, y, z transformiert und erst 
jetzt geglattet. Dieses Verfahren erfordert etwas mehr Rechenmittel als die 
bisher genannten, aber es vereinigt im tibrigen die Vorteile von 4a und 4c, 
ohne deren Nachteile aufzuweisen: Die Differentiation kann unmittelbar an ~ 
den Visierwerten vorgenommen werden, und die Glattung von v erfolgt in 
einem raumfesten System. Beztiglich der erreichbaren Genauigkeit ist dieses 
Verfahren als das beste zu bezeichnen. 


e) Andere Darstellungen von v 


Viele Gerate verzichten darauf, die Geschwindigkeit in einem eigentlichen 
Koordinatensystem darzustellen, und bilden lediglich aus den Visierwerten die 
Ableitungen «, A, ¢. Die Addition zum Messpunktvektor erfolgt dann mit Hilfe 
verschiedener zum Teil approximativer Verfahren, die im nachsten Abschnitt 
geschildert werden. Solche Lésungen erméglichen in den meisten Fallen eine — 
bedeutsame Materialersparnis, sind aber im Sinne der Glattung ungiinstig. 


6. Addition M+ ut 


Die Multiplikation des Geschwindigkeitsvektors v mit der Flugzeit t ergibt — 
den Vorhaltevektor. Dieser Prozess wird mit normalen Rechenelementen ver- 
wirklicht. Er bildet fiir das Studium der Koordinatensysteme nichts Interes- 
santes und wird daher hier nicht beschrieben. Der entstehende Treffpunkt- — 
vektor muss in ein Koordinatensystem transformiert werden, welches fiir die 
Bildung der ballistischen Funktionen geeignet ist (vgl. Figur 3). Meistens sind 
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dies Zylinderkoordinaten («;, e,7, hp) oder Kugelkoordinaten (a7, Ar, er). Nun 
tberstreichen die Komponenten von Mund von Teinen sehr grossen Bereich. 
Im Beispiel einer 7,5-cm-Batterie liegen die kartesischen Messpunktkoordi- 
naten x und y etwa im Bereich von +. 15.km und h zwischen 0 und 8 km. Dies 
stellt an die Prazision der Koordinatentransformatoren grosse Anforderungen, 
und bei gegebenen relativen Fehlern der Rechenelemente verursacht eine 
solche Transformation einen Genauigkeitsverlust fiir das ganze Gerit. Es ist 
wesentlich besser, wenn gemiss Figur 4 der Messpunktvektor zwischen Visier 
und Ballistik ohne Transformation iitbertragen wird. Die Uberfiihrung von M 
in T erfolgt durch Addition der Differenzen Ax, 4A, Ae; diese Addition wird 
mit Hilfe von Differentialgetrieben ausgefiihrt und verursacht im Gegensatz 
zu allen andern Rechenoperationen keinen Genauigkeitsverlust, da in solchen 
Getrieben der ganze Bereich durch beliebig viele Umdrehungen dargestellt 
werden kann. Die durch die iibrigen Rechenelemente verursachten relativen 
Fehler sind nur noch in Aa, AA, Ae enthalten und daher wesentlich kleiner als 
in Figur 3. 

Nicht in allen Geraten verwendet man die vorteilhafte Anordnung von 
Figur 4; es gibt auch Zwischenlésungen, welche von den drei Visierkoordinaten 
nur « direkt zur Ballistik weiterleiten, die ibrigen jedoch transformieren. 


a) Addition in kartesischen Koordinaten 


Wenn sowohl Messpunktvektor als auch Geschwindigkeitsvektor in orts- 
festen kartesischen Koordinaten gegeben sind, so erfolgt die Bildung von T 
in einfachster Weise nach folgenden Gleichungen: 


tp = HVT, Vr H+ OTS hp Shay, 7. 


Dieses Verfahren ist iibersichtlich und erfordert wenig Rechenelemente, 
bringt aber den oben dargelegten bedeutenden Genauigkeitsverlust mit sich 
und wird daher heute kaum mehr verwendet. 


b) Gemischte Addition 


Haufig werden die aus dem Visier kommenden Messpunktkoordinaten «, 
4, e in Zylinderkoordinaten («, e,, 1) transformiert. Diese Transformation be- 
trifft nur A und e, wahrend « unverandert bleibt. Zu diesen Koordinaten werden 
die aus der Geschwindigkeit gewonnenen Differenzen Aa, Ae,, Ah addiert. Dies 
stellt eine Zwischenlésung zwischen Figur 3 und 4 dar. Genauigkeitsmassig ist 
das gewonnene Resultat besser als in Figur 3, da e, und h einen kleineren 
Bereich iiberstreichen als x, y, . Der verbleibende Fehler ist aber immer noch 
erheblich grésser als im Verfahren nach Figur 4. — Fiir die Bildung der drei 
Differenzen werden drei verschiedene Verfahren verwendet, wobei der Ge- 
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schwindigkeitsvektor gemass Abschnitt 5b gewonnen wird. Fiir da und Ae, 
wird im Gerat eine geometrische Abbildung gemass Figur 1 ausgefiihrt. Aus 
dem dargestellten Dreieck wird Aa direkt abgenommen und durch ein Diffe- 
rentialgetriebe zu « addiert. Ae, erscheint als eyr — & und wird meist mit 
Hilfe eines Stahlbandes gebildet. 4h entsteht in der Form v, t und wird in ge- 
eigneter Weise zu h addiert. Es ist ersichtlich, dass dieses Verfahren nur in 
bezug auf « den Vorteil der erhdhten Genauigkeit aufweist. 


c) Addition gemdss Figur 4 


Falls v entsprechend Abschnitt 5c in einem bewegten kartesischen System 
vorliegt, so kénnen daraus mit Hilfe einfacher Operationen die Werte Aa, AA, | 
Ae gebildet und zu den Messpunktkoordinaten addiert werden. Dadurch sind 
alle Méglichkeiten zur Reduktion des Einflusses, den der relative Fehler der 
Rechenelemente verursacht, ausgeniitzt. Wird im Interesse einer einwand- 
freien Glattung der Geschwindigkeitsvektor nach Abschnitt 5d in einem raum- 
festen kartesischen System gebildet, so konnen aus den Komponenten v,, vy, v;, 
ebenfalls die Differenzen Ma, 4A, Ae dargestellt werden, immerhin mit einem | 
erheblich vergrésserten Rechenaufwand. Dieses System vereinigt als einziges 
die Vorziige, die durch die Anordnung von Figur 4 gegeben sind, mit denjenigen 
der konstanten Geschwindigkeitskomponenten. 


d) Direkte Bildung von Aa, Ad, Ae 


Durch geeignete trigonometrische Formeln lassen sich die Differenzen 
Aa, Ad, Ae direkt aus den Visierkoordinaten «, A, e und ihren zeitlichen Ablei-— 
tungen bilden. Dieses Verfahren ist allgemein recht giinstig. Immerhin ist zu 
beachten, dass v als Vektor tiberhaupt nicht in Erscheinung tritt. Die Tat- 
sache, dass diese Grdsse durch den Feuerleitenden nicht abgelesen werden : 
kann, ist als bedeutender Nachteil zu betrachten. Damit sind natiirlich auch 
die weiter oben geschilderten Einschrankungen beziiglich der Glattung ver- 
bunden. | 


e) Taylor-Entwicklungen 


Oftmals besteht der Wunsch nach einem Feuerleitgerat, das mit einem 
minimalen Aufwand an Rechenelementen brauchbare Resultate liefert. Hier 
greift man oft zur Taylor-Entwicklung und stellt die Differenzen wie folgt dar: 

: 


Aa=at+R,, AM=Ar+ Ry, de= ert Ry. 


R,, Ry, Rz sind die Restglieder. Ihr Anteil an der gesamten Summe wird dann 
besonders gross, wenn die Visierwerte eine hohe zweite Ableitung besitzen, 
was in der Umgebung des Wechselpunktes der Fall ist. Die approximative 
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Darstellung dieser Restglieder durch méglichst einfache trigonometrische For- 
meln, welche nur die Visierwerte und ihre ersten Ableitungen enthalten, bildet 
fiir den Mathematiker eine dankbare Aufgabe, und es sind zahlreiche verschie- 
dene Systeme aufgestellt worden, deren-Genialitat zum Teil bewundert werden 
muss. Alle diese Systeme geben nur angendherte Werte, deren Fehler aber, 
mit Ausnahme extremer Fille, gegeniiber der ballistischen Streuung in den 
Hintergrund treten. 
Es ist auch versucht worden, die Restglieder durch Zufiigen einer zweiten 
und dritten Ableitung darzustellen, also zum Beispiel 


a ee ae 

Diese Methode scheitert aber daran, dass bei mehr als einmaliger Differentiation 
der Visierwerte die immer vorhandenen Stérungen unzulassig verstarkt wer- 
den. Dies fallt bei grossen Flugzeiten um so mehr ins Gewicht, als die héheren 
Ableitungen noch mit héheren Potenzen der Flugzeit zu multiplizieren sind. 
Eine fiir Kleinkaliberwaffen allenfalls noch gangbare Darstellung ist die fol- 
gende: : . 
Aa=att+ate 

(Das Restglied enthalt die erste Potenz der Flugzeit und nicht die zweite, wie 
es fiir die Taylor-Entwicklung richtig ware). Durch geeignete Wahl der Kon- 
stanten c kénnen die Fehler, welche ohne Restglied im Wechselpunkt auf- 
treten wiirden, erheblich vermindert werden. 


f) Richtgerdtesteuerung 


Allgemein ist zu sagen, dass die vollstandige Taylor-Entwicklung nicht die 
erstrebenswerte Arbeitsweise fiir ein Feuerleitgerat darstellt, selbst wenn man 
von Stérungen absieht. Der Satz, dass eine Funktion vollstandig bestimmt ist, 
wenn an einem einzigen Punkt Funktionswert und alle Ableitungen gegeben 
sind, gilt nur fiir analytische Funktionen, zu denen jedenfalls der Weg eines 
Flugzeuges nicht gehért. Um iiber den zukiinftigen Flugweg mdglichst genaue 
Angaben zu erhalten, ist es vielmehr nétig, méglichst weit in die Vergangenheit 
zuriickgreifend die Vorgeschichte zu erfassen. Diese fiihrt zu einer Integral- 
darstellung des Flugweges und damit zu Geraten, die vom geschilderten ‘Autf- 
bau erheblich abweichen. 


7. Darstellung der Ballistik 


Dem Geschiitz miissen die drei Werte «7, ¢, t zugefiihrt werden. Der Treff- 
-punkt T muss daher zunachst in Zylinderkoordinaten («;, ¢7, 7) oder Kugel- 
koordinaten (a7, Ar, er) transformiert werden. Die Treffpunktseite a; kann 
dann direkt, ohne ballistische Umrechnung, dem Geschiitz zugefiihrt werden 
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(von der Derivation, welche in jedem Fall nur eine kleine Korrektur verur-_ 


sacht, wird hier abgesehen). Die verbleibenden zwei Komponenten liegen in 
der Schussebene; im Fall e,7, 4p handelt es sich um kartesische Koordinaten in 


dieser Ebene, im Fall 4;, er um Polarkoordinaten. Diese Komponenten miissen ~ 


durch den ballistischen Auswerter in die Grdssen ¢,t umgerechnet werden. 
Je nach der Konstruktion des ballistischen Auswerters eignet sich das eine 
oder das andere System besser. Daraus ist ersichtlich, dass die Natur des 
ballistischen Rechners auf die Wahl aller vorangehenden Koordinatensysteme 
einen Einfluss haben kann. 


a) Einfache ballistische Korper 


Jede der ballistischen Gréssen ist eine Funktion zweier Variablen, also etwa 


im Fall der Koordinaten e,7, hr: 


E= Ep, hp), T= T(&z, Az) - 


Es ist naheliegend, diese Funktion mit Hilfe von zwe: Kurvenkorpern zu er- | 
zeugen. Dies sind graphische Darstellungen der Funktion von zwei Verander- — 


lichen in Form einer Héhenkarte, wobei ein Fiihler in zwei Richtungen bewegt 
werden kann und den gesuchten Funktionswert als Héhe abtastet. Mittels 


solcher Kurvenkérper kann — mit gewissen Einschrankungen — jede Funktion — 
erzeugt werden, und die Fabrikation benétigt keine mathematische Vorarbeit. — 
Die Korper sind aber tiberaus empfindlich und teuer und erméglichen nur eine ~ 
beschrankte Genauigkeit ; dabei sind in der erwahnten Darstellung die Genauig- — 


keitsanforderungen sehr hoch, geht doch der Bereich zum Beispiel fiir ¢ von 


0-1600 A°/o), wobei ein Fehler von nur wenigen A°/,, zugelassen werden kann. © 


b) Abspaltung linearer Summanden 


Fiir das Verfahren, das im folgenden beschrieben werden soll, eignet sich | 


die Darstellung des Treffpunktvektors durch Polarkoordinaten (A;, er) in der 


Schussebene besonders gut. Es sind dann folgende Funktionen zu bilden: 
&=e(Ap, ep), t= T(Az, er). 
Es ist sofort ersichtlich, dass « vorwiegend von A; abhangt, wahrend 7 in erster 
Linie durch die Entfernung e; bestimmt ist. Die Funktionen werden daher in 
der folgenden Art aufgespalten: 
é=A+o(d;,e7), t= om + flap, ep) . 


Die Schusswinkelfunktion o(A;, e-) durchlauft nun einen wesentlich kleineren 


Bereich und muss daher mit kleinerer relativer Genauigkeit erzeugt werden 
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was besonders bei Verwendung von ballistischen K6rpern ein Vorteil ist. v,, 
ist eine Konstante, welche der mittleren Geschossgeschwindigkeit entspricht. 
Sie ist so gewahlt, dass die Funktion /(A;, e7) méglichst klein wird und daher 
ebenfalls mit geringer relativer Genauigkeit erzeugt werden kann. (Die Multi- 
plikation der Entfernung e; mit der Konstanten 1/v,, geschieht durch ein 
Getriebe und verursacht keinen Genauigkeitsverlust.) 

Das beschriebene Verfahren der Zusammensetzung der ballistischen Werte 
aus einem grossen, exakt darstellbaren Anteil und einer kleineren, nur mit 
beschrankter Genauigkeit zu bildenden Korrektur bietet beziiglich Aufwand 
und Kosten eines Gerates bedeutende Vorteile. 


c) Zuriickfiihrung auf Funktionen einer Variablen 


Die Erzeugung von Funktionen zweier Variablen, wie o(Ar, er) oder f(A, e7) 
im vorigen Abschnitt, erfordert Kurvenkérper. Solche Kérper sind mit Nach- 
teilen behaftet: Sie sind teuer in der Fabrikation, beanspruchen viel Raum 
und sind empfindlich gegen Erschiitterungen. Demgegeniiber gibt es mehrere 
Vorrichtungen zur Erzeugung von Funktionen mit einer Variablen, die in Her- 
stellung und Verwendung iiberaus vorteilhaft sind; es sind hier zu nennen: 
Kurvenscheiben, Zugmittelgetriebe, Gelenkwerke. Besonders die letzteren, die 
nur Kurbeln und Gelenke enthalten, sind sehr genau und robust [3]. In elektro- 
mechanischen Geraten kénnen dazu noch Funktionspotentiometer und Funk- 
tionsdrehkondensatoren verwendet werden. 

Daher besteht der Wunsch, die gesamte Ballistik nur mit Funktionen einer 
Variablen darzustellen, welche mittels Additionen und Multiplikationen zusam- 
mengesetzt werden. Beispielsweise ist folgender Ansatz mdéglich: 


(Ay, er) = A, (A) E,(er) + Ag(Ar) + Eo(er7) - 


Der Entscheid, ob dieser Ansatz geniigt oder eventuell sogar vereinfacht wer- 
den kann, und das Auffinden der vier neuen Funktionen bilden eine umfang- 
reiche mathematische Arbeit, wobei zweckmassig mit den Methoden der Sta- 
tistik vorgegangen wird [4]. Der Erfolg einer solchen Darstellung hangt in ent- 
scheidendem Masse davon ab, ob die Treffpunktwerte in einem zweckmassigen 
Koordinatensystem gegeben sind. Die Wahl der Koordinatensysteme in einem 
Feuerleitgerat ist daher eng mit der Art der Darstellung der Ballistik verkniipft. 

Es ist darauf hinzuweisen, dass eine solche Darstellung die vorgegebene 
Funktion nur approximativ erzeugen kann. Der Fehler lasst sich aber durch 
geeignete Wahl des Ansatzes und durch Hinzunehmen gentigend vieler Glieder 
beliebig klein machen, so dass er schliesslich gegentiber den Ungenauigkeiten 
der Rechenelemente in den Hintergrund tritt. itis 

Falls die gewahlten Rechenelemente fiir die Erzeugung von beliebigen, vor- 
gegebenen Funktionen ungeeignet sind und nur Additionen, Multiplikationen 
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und Divisionen ausfiihren kénnen — dies ist zum Beispiel bei der Wahl eines 
digitalen Rechenwerks der Fall —, so kénnen die Funktionen etwa durch Poly- 
nome oder durch gebrochene rationale Funktionen approximiert werden. Die 
Erfahrung zeigt, dass sich ganz besonders die letzteren fiir die Darstellung bal- 
listischer Beziehungen gut eignen. 


8. Korrekturen 


a) Tagesunstimmigkerten 


Durch Abweichung des Luftgewichtes 6 und der Geschossanfangsgeschwin- 
digkeit v, von den Wertén, die bei der Errechnung der Ballistik zugrunde 
gelegt wurden, entstehen Anderungen in der Flugbahn. Der Wind wird in 
dieser Betrachtung ausser acht gelassen. Die Natur der nétig werdenden Kor- 
rekturen, mit denen ¢ und t versehen werden muss, ist so kompliziert, dass in | 
jedem Fall nur approximative Lésungen méglich sind. Eine haufig verwendete 
und iibersichtliche Darstellung ist die folgende: | 

Ae=f,(T) A6 + f,(T) dv, Ar =f,(T) Ad-+ f(T) Avy 

Die Funktionen /,, ... , /, sind skalare Funktionen des Vektors T, wobei T 
in Zylinder- oder Kugelkoordinaten dargestellt ist ; die Funktionen sind natiir-' 
lich unabhangig von «; und haben daher nur zwei unabhangige Veranderliche. 
In diesen Formeln treten die Tagesunstimmigkeiten 46 und Av, nur linear auf. 
Dadurch ist insbesondere der gegenseitigen Beeinflussung nicht Rechnung ge- 
tragen, und es treten bei grossen Tagesunstimmigkeiten erhebliche Fehler auf. 
Trotz dieser Unzulanglichkeiten bedingt das beschriebene Rechenschema einen 
bedeutenden Materialaufwand, indem 4 Kurvenkérper (oder analoge Vorrich- 
tungen) und 4 Multiplikatoren nétig sind. 

Ein wesentlich eleganteres Verfahren erhalt man durch Anbringen der Kor- 
rekturen vor dem ballistischen Rechner, also nicht an den ballistischen Gréssen, 
sondern an T. Dies entspricht der Bildung einer fiktiven Treffpunktverlagerung 
AT nach folgender Beziehung: 


AT = S,(T) 46 + S,(T) Av,. (3) 


T und AT haben auch hier nur zwei wesentliche Komponenten, da das Azimut 
ar keine Rolle spielt. S, und S, sind vektorielle Funktionen des Vektors T. 
Dem ballistischen Auswerter wird der Vektor T+ AT zugefiihrt. Dieses Ver- 
fahren ist, ebenso wie das vorher geschilderte, nur approximativ. Es hat aber 
den Vorteil, dass die Tagesunstimmigkeiten in ¢ und t nicht linear zur Auswir 
kung kommen, was eine bessere Anndherung an die Wirklichkeit bedeutet. Be- 
sondere Vereinfachungen im Gerat bringt die folgende Modifikation: Der 
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/ektor AT hat nur eine e-Komponente und hangt nur von e ab, so dass S, 
ind S, skalare Funktionen einer Variablen werden..Um den dadurch entste- 
ienden Fehler auszugleichen, rechnet die Ballistik wie folgt: 


Fir die Bildung von eT) = T+ yu, AT. | 


4 
Fir die Bildung vont: T,= T+ yu, AT. | e 


t, und jl, sind dimensionslose Zahlen. Eine davon kann willkiirlich gleich Eins 
esetzt werden, und die andere ist durch Versuch zu bestimmen. Das System 
4) braucht nur noch zwei Funktionen einer Variablen, gegeniiber vier Funk- 
ionen zweier Variablen im Falle (3). 


b) Ladeverzug 


Der Ladeverzug k entsteht dadurch, dass zwischen dem Tempieren und 
lem Abfeuern des Geschosses eine gewisse Zeit — etwa 2-3 s — verstreicht. Die 
fempierung eines Geschosses, das zu einer bestimmten Zeit tempiert wird, ist 
laher nicht gleich wie die Flugzeit eines Geschosses, das zur gleichen Zeit abge- 
euert wird. Die getrennte Berechnung von t und # erfordert, wenn sie exakt 
tfolgen soll, im Gerat einen erheblichen Mehraufwand, indem die Addition 
les Vorhaltevektors zweimal erfolgen muss. Auf Kosten gewisser Ungenauig- 
‘eiten kénnen aber bedeutende Einsparungen erzielt werden. Wohl die ein- 
achste Berechnung der Tempierung # aus der Flugzeit erfolgt nach folgender 
Seziehung: 

B=t+Ttk. 


Jiese Formel ist hauptsadchlich dann ungenau, wenn Tt eine hohe zweite Ablei- 
ung besitzt. 


c) Parallaxe 


Die Parallaxe kommt dadurch zustande, dass Feuerleitgerat und Geschiitz 
licht am gleichen Ort stehen; dies bedeutet, dass der errechnete Treffpunkt 
m Sinn einer Translation in ein neues Koordinatensystem umgerechnet werden 
nuss. Die Addition des Parallaxvektors P zum Treffpunktvektor T ist véllig 
nalog der Addition M+ ut, und in der Regel wird in einem Gerat auch fir 
yeide Aufgaben das gleiche Verfahren angewendet. Es kann daher fiir die ver- 
chiedenen Méglichkeiten auf die Abschnitte 6a, 6b und 6c verwiesen werden. 
a der Parallaxvektor in der Regel kleinere Werte annimmt als der Vorhalte- 
ektor, ist es hier méglich, in vermehrtem Masse von Naherungslésungen 
tebrauch zu machen. Besonders im Verfahren nach 6b gelingt es, fiir Aa, AA, 
le sehr einfache, wenn auch nur approximative Ausdriicke zu finden. Ein 
enaues und trotzdem sparsames Verfahren lasst sich anwenden, wenn gemass 
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Abschnitt 6 c der Vorhaltevektor v t in kartesischen Koordinaten dargestellt 
ist. Dann wird der Patallaxvektor direkt komponentenweise dazu addiert, und 
die Umrechnung auf Aa, AA, Ae erfolgt fiir beide gemeinsam. 


9. Statistischer Aspekt 


An verschiedenen Stellen dieser Arbeit wurde erwahnt, dass fiir gewisse 
komplizierte Umformungen Naherungslésungen gefunden werden kénnen. Die 
Verwendung solcher Lésungen fiihrt dazu, dass ein Gerat nicht nur zufallige 
Fehler aufweist, die durch ungenaue Herstellung und Justierung verursacht 
werden, sondern dass auch noch systematische Fehler entstehen, also solche, 
die fiir den gleichen Gefechtsablauf immer und bei jedem Exemplar des Ge- 
rites gleich sind. Es ist denkbar, dass es Gefechtsablaufe gibt, in denen sich 
zeitweise die Ungenauigkeiten verschiedener Rechenelemente summieren, so 
dass in diesen Augenblicken ein Schussfehler entsteht, der so gross ist, dass 
die Trefferwahrscheinlichkeit fast auf Null reduziert wird. Es ware aber falsch, 
daraus die Unbrauchbarkeit des betreffenden Gerates ableiten zu wollen. Um 
ein Mass fiir die Genauigkeit eines Feuerleitgerates zu erhalten, muss vielmehr 
eine mittlere Treffererwartung tiber alle moglichen Flugwege gebildet werden, 
wobei die verschiedenen Flugwege je nach ihrer Haufigkeit mit verschiedenen 
Gewichten zu versehen sind. Es ist zu betonen, dass oft gerade die mit gewissen 
systematischen Fehlern behafteten Rechenverfahren auf Lésungen fiihren, die 
infolge der Einfachheit und Zuverlassigkeit der Rechenelemente in der Mehr- 
zahl der Gefechtsablaufe besonders genaue Resultate liefern. 
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Summary 


The best choice of co-ordinate systems in fire control computers depends not 
only upon the type of computing elements used, but also upon the formulae b : 
which the geometric and ballistic computations are to be done. By cert 
selecting co-ordinate systems at every point of the computing process, arnelberd : 


simplifications may be made which will increase both iabili 
“enna accuracy and reliability 


(Eingegangen: 12. Dezember 1954.) 
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Das Vertahren der abgekiirzten Iteration fiir algebraische 
Eigenwertprobleme, insbesondere zur Nullstellenbestimmung 
eines Polynoms 


Von FRIEDRICH L, BavEer, Miinchen?) 


Das klassische Iterationsverfahren von DANIEL BERNOULLI [5]*) zur Be- 
stimmung der absolut gréssten Nullstelle eines Polynoms kann bekanntlich 
als Vektoriteration [11], durchgefiihrt mit der Frobeniusschen begleitenden 
Matrix des Polynoms, aufgefasst werden (AITKEN [1], Fry [7]). Es wird in dieser 
Arbeit gezeigt, dass die Berechnung der (linear konvergenten) Folge der Itera- 
tionsvektoren so abgekiirzt werden kann, dass eine quadratische Konvergenz 
erreicht wird. Dabei bleibt dieses im folgenden « abgekiirzte Iteration » genannte 
Verfahren gegen Rundungsfehler selbstkorrigierend. 

Die Tragweite der abgekiirzten Iteration zeigt sich jedoch besonders bei der 
Bestimmung weiterer Nullstellen. Die gewéhnliche Vektoriteration an einem 
passend orthogonalisierten Ausgangsvektor ist in diesem Fall bekanntlich 
numerisch instabil. Die abgekiirzte Iteration jedoch hat gleichzeitig stabilisie- 
renden Charakter, es ist also keinerlei laufende Orthogonalisierung no6tig. Sie 
erfordert keine Abanderung des Ausgangspolynoms, die Verschleppung von 
Rundungsfehlern ist also verhindert und eine selbstkorrigierende Berechnung 
aller Nullstellen erméglicht. Gegen Rechenirrtiimer schiitzt eine einfache Ver- 
probungsméglichkeit. 

Die algebraischen Grundlagen des Verfahrens wurden in einer gesonderten 
Abhandlung vorweggenommen [4]. Der Schwerpunkt liegt in der vorliegenden 
Arbeit auf den Gesichtspunkten, die fiir die praktische Durchfiihrung gelten; 
jedoch ist alles fiir das Verstaéndnis Wesentliche aus einem vorangestellten 
Abriss der Theorie zu entnehmen. Wir behandeln in dieser Arbeit auch die der 
besseren Einsicht in das Verfahren dienende Verallgemeinerung auf Matrizen- 
Eigenwertprobleme. Fiir praktische Zwecke wird allerdings oft, insbesondere 
wenn Eigenvektoren nicht gesucht sind, ein direktes Verfahren, das unmittel- 
_ bar auf das Sikularpolynom zielt, weniger Aufwand erfordern. 

Das hier geschilderte Verfahren wurde in umfangreichen Rechnungen nu- 
merisch erprobt. Ich darf Herrn Prof. A. WALTHER, Darmstadt, und der Deut- 
 schen Forschungsgemeinschaft, die dies ermoglichten, herzlich danken. 


1) Mathematisches Institut und Arbeitsgruppe fiir elektronische Rechenanlagen der Tech- 


nischen Hochschule Miinchen. 1 hee 
2) Die Zahlen in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 32. 
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§ 1. Bezeichnungen 
Wir bezeichnen mit 


P(x) Biss Re Ee (1) 


das Polynom, dessen Nullstellen bestimmt werden sollen. Die Koeffizienten 
a, mogen beliebig, reell oder komplex sein. Es kann jedoch fiir numerische 
Aufgaben stets angenommen werden, dass a) = 1_und 4, + O ist. Die Null- 
stellen seien nach ihren Betragen numeriert, 


(| Sle ele = oe (2) 


Der Einfachheit halber sei zunaichst angenommen, dass alle Nullstellen einfach 
sind. 
Fiir ein beliebiges Polynom Q(x) gibt es dann die Partialbruchzerlegung 


Aden oa LCI PS EL a 23 eas 


wo H(x) eine ganze Funktion ist. 
Summen sollen, wo nichts anderes gesagt ist, stets von 1 bis  laufen. Mit 


Pa (a) = eat 


m=1,...,” Come 


a ef P(x 
Pr) = ee (5) 
die normierten. Es ist wegen P’(x) = X P,(x) 
Pe) ie aca) (6) 
und also 
Q(x) = 37 OG) P,(x) mod P(x) . (7) 


Kongruenzen sollen in dieser Arbeit stets im Restklassenring mod P(x) ver- 
standen werden. 


Nach der Definition von P,,(x) ist 


SP, (acres ee 
und also i‘ a y 


# O(n) = SE OE) Pia. (9) 


Das Polynom vom Grade n — 1 auf der rechten Seite sei als Q(x) bezeichnet, 


Q(x) SE O(E) B(x) . (10) 
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Der Koeffizient von x"~1 in Q‘(x) ist 


eine gewichtete Potenzsumme, Mit Q(x) sei das normierte Polynom 


ee Q(x) 


oe 


OM(x 


worin der Koeffizient von x"~! gleich 1 ist, bezeichnet. 


(12) 


§ 2. Theoretische Grundlagen des Verfahrens der abgekiirzten 
Iteration 


1. Die Iteration mit x mod P(x) 


Die Polynome Q(x) oder Q(x) kénnen rekursiv berechnet werden, 


QD (x) = x Q(x) — a, P(x) , (13) 
OFN(x) = Le ons) — P(x) ; (14) 


Wir werden diesen Prozess kurz « (gewéhnliche) Iteration mit x mod P(x) » 
nennen. 
Falls 
BE) = QE s=* OC a) = O4E,) + Ors, | (15) 
und 


|é,| a | Epcy 


iiberwiegt der Einfluss von & in Q" den aller tibrigen Nullstellen, wenn 7 unbe- 
schrankt wachst; man hat 


Tn, OF (x) PAE (16) 


1—>0o 


Die Iteration mit x mod P(x) fiihrt also dann direkt zum Polynom P,(x), 
von dem die Nullstelle &, abgetrennt ist, und gleichzeitig auch zu &, selbst: 
E, ist der Faktor, mit dem sich Q")(x) in der Grenze multipliziert, 


Q(x) —£, OP) > 0. (17) 


Die Konvergenz ist linear, der Fehler wird im Limes bei jedem Schritt mit 


&.\/|&..,]| multipliziert. : 
| ee Prozess liegt implizit die von LAGRANGE [10] herriithrende Variante 


der Bernoullischen Methode*) zugrunde. 


3) Siehe [4], § 1.2. 
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Die Iterationspolynome Q(x) spielen bei LAGRANGE nur eine Nebenrolle. 
Gerade sie sind es jedoch, die eine betrachtliche Abkiirzung der Iterationsfolge 


gestatten. 


2. Abkiirzung der Iteration 


Auf die Frage, ob man Q(x) unmittelbar aus Q(x) erhalten kénne, lautet 


namlich die Antwort: 
Es ist: 


QP) = YL OE) Pix), (WMP =D (E 0G)" 2 » 


da 
B(x) Px) = P(E) P(x) = 6, PQ) j 
also 
0° (x) =[Q°W)F 
falls 


OG) = 0 oder 1 (vy =) we a. (18) 


Zur expliziten Berechnung von Q©')(x) bieten sich verschiedene Méglich- 
keiten, um [Q(x)]?, ein Polynom vom Grad 2 — 2, modulo P(x) zu reduzieren. 
Am iibersichtlichsten ist die folgende*). Sei 


OMe) = diye (19) 

dann ist 
[Ow] = Sy 2 Oe) = Fy O(a). (20) 
Es sind also durch gewoéhnliche Iteration die Polynome Q°*”,..., Q@*"-) 


herzustellen, dann erfolgt eine gewichtete Summierung. Die Gewichtsfaktoren 
liefert das Polynom Q(x). 

Da bei jedem dieser Schritte der Iterationsindex verdoppelt wird, ist der 
Fehler nach jedem Schritt von der Ordnung des Quadrates des vorherigen 
Fehlers, die Konvergenz ist nunmehr quadratisch. 


3. Wahl des Ausgangspolynoms 


Zur Bestimmung von é,, falls |&,| > |é,., 


, hat man also die Bedingungen 
(15) und (18) zu erfiillen. Dazu geniigt es, 


Q(x) =79(y “1 — SP = SPW (21) 
zu wahlen. é 


Der erste Lauf der Iteration ist insbesondere stets mit 10(x) = 1 zu be- 
ginnen. 


4) Weitere siehe [4], § 2.1 und § 2.3. 


' 
: 
: 
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Wird sodann beim r-ten Lauf der abgektirzien Iteration, startend mit 
Q(x), P,(x) bestimmt, so kann der neue, (ry + 1)-te Lauf mit 


r+10)(x) * "O(x) =. P (x) (22) 


unmittelbar angeschlossen werden. Dies geht so lange gut, als keine Nullstellen 
gleichen Betrags an die Reihe kommen. 


§ 3. Nullstellen gleichen Betrags 
1. Die Bildung der abgeleiteten Iterationsfolge 


Wir beginnen mit dem Fall, dass zwei Nullstellen von gleichem Betrag sind. 
Dies trifft bei reellen Polynomen zu, wenn es sich um ein einfaches konjugiert- 
komplexes Paar handelt. 

Ausgangspolynom sei Q(x), es sei |&,| = |&4,| > |&45|. Die Folge der 
Iterationspolynome ’Q)(x) konvergiert nicht. Ist &, = &,,, = @ e’”, so oszilliert 
sie in allen Gliedern. Jedoch konvergiert nun die Folge der « abgeleiteten » 
Polynome 30“ (x) vom Grad m — 2, die man durch einen Prozess erhilt, der 
formal iibereinstimmt mit dem der « cross-multiplication» von RovurtH [12]. 

Es sei 


(23) 


5Q(x) sei wieder so normiert, dass der Koeffizient von «”~? Eins ist, man hat 


also Lage a 
7O"%F D(x) ss 70 (x) 


rol (x) = 24 
20" (2) Koeffizient von "~~? a 
Es gilt sodann?®) 

lim £0 (x) = P, ,41() (25) 


wo P, ,,, das Polynom vom Grad m — 2 ist, von dem der zu é, und é,,, gehorige 
quadratische Faktor 


Sy e43\) = (x - Col (x = S54) = x* — p x+q 
_ abgetrennt ist: 
: ee ee ; (26) 


ih) 


Damit ist auch S, ,..;(x) bestimmt. Im tibrigen gilt im Limes 


; re Nees ae p LY gt ae q 79") > On (27) 
ye 5) Spezialfall aus [4], § 3.2. Der dortige Beweis tibertragt sich sofort, wenn man bedenkt, dass 
Ei: 7Q(x) mit P(x) die Linearfaktoren x — Eye CeO — &,_4 gemeinsam hat. 
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In dem bisher ausgeschlossenen Fall mehrfacher Nullstellen, hier ESeray 
konvergieren die 3Q(x) ebenfalls®) linear, und zwar (vgl. Abschn. 3) erheblich 
rascher als die logarithmisch konvergenten ’Q")(x). 


2. Gewinnung des neuen Ausgangspolynoms 


Der Start eines neuen Laufes nach der Bestimmung von é, und &,,, erfor- 
dert die Bildung von *+2Q(x) = Q(x) — [P.(x) + P.,,(x)]. Das zu subtrahie- 
rende P(x) + P.,,(x) kann fiir reelle Ausgangspolynome unter Vermeidung 
einer Rechnung im Komplexen so bestimmt werden: Es ist 


rR HN a & =. é, 
Pa) + Preals) = [=peG® + poe ty | Porn) (28) 


(Fy) 


Der zweite Faktor hat sich durch die Iteration direkt ergeben. Ausserdem 
ist &.+ &,,, = p und &,&,,, = q bekannt. Mit ihrer Hilfe lasst sich die eckige 
Klammer umformen, es ergibt sich 


[~— (p/2)) +H 


P,(x) + Pao) a [q — (p2/4)] G2? + H? Py rai (%) Ne (29) 


WO 
(cee False) — Fr rsalFrsa) ; TT PaaS) Tr Pet een (30) 
=(é, ~ eeu) 2 
ist. Fihren wir die Koeffizienten m, von P,,,, ein, 
Pn 41(%) == lg hele ate yaa nae + Ty, 9 , (31) 
so wird 
N—2 1 n—2 
G= 2 Ny 1 Uy, 9» » H= Tos (ine J Ny 2) My, 2» > (32) 
Dabei ist 
Ny = ae - uae 
r Sr+1 


eine totalsymmetrische Funktion y-ten Grades, rekursiv definiert durch 


1 —=— pve eel al eee ; (33) 


= 1 
NM=1, 7 1=0, Dao ae (34) 


Die angegebene Form hat den Vorzug, dass sie im G 
lich bleibt, es geht 6 ie im Grenzfall &, = &,,, end- 


Co HFG) de Pa ley (35) 


*) Vgl. [4], § 4. Die dortige Bedingung e{) + 0 ist sinngemass erfiillt 
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Im ibrigen sind mannigfache Umformungen méglich. Mit Hilfe von (33) erhalt 
‘man zum Beispiel 


P,(x) a +1(%) 


he TO os “— be v7 can 
(5 iain (5 "sya (36) 


j 
oo / 


vias n—2 7 m2 ee °| Fr re1(%) « 
| ai(& ny m-1-+] 85 Ny_9 sh ks| Per (5 Ny Tot) | 
0 v0 / 0 


3. Bildung wetterer abgeleiteter Iterationsfolgen 


Im Falle, dass m > 2 Nullstellen vom gleichen Betrag, jedoch sdmtliche 
verschieden sind, 


irs |é,| a, [Enea| ie Ee oor er a Spit , 


konvergiert erst die m-te der abgeleiteten Folgen 70 (x) von Polynomen vom 
Grade n — x, die rekursiv durch fortgesetzte Durchfiihrung der « cross-multi- 
plication » gebildet sind’) 


dp rov+) rol) 
5" (x) ee es x12 — (ee Borie it) x (37) 
Koeffizient von x” * 
Es gilt’) 
‘ sal P(x) 
ben thrial a De ' = : 38 
jim,,Q (x) r,rt+l,...,r+m—1 Se ha: fey OA) ( ) 


es entsteht auch hier direkt das Polynom, von dem der Faktor 
S(x) ia yada Es aa (x -; é,) (x m7 ft) er (x ee (39) 


— dessen Nullstellen die gesuchten sind — abgetrennt ist. Deren Bestimmung 
bietet nun bekanntlich keine Schwierigkeiten mehr, mit Hilfe der gebrochen 


linearen Substitution ’ 
a+ a 


oe 


‘die den Kreis vom Radius 9 auf die reelle Achse abbildet, entsteht aus S(x) ein 


Polynom 
; S(z) = (2 — 4)" S(o 25) 


he, 
Bit reellen Nullstellen, das fiir reelles S(x) entweder gerade oder ungerade ist. 


7) Vel. [4], § 3.2. 
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Zur Bildung des neuen Ausgangspolynoms 
EOL E > Q(x) ma [ P,(x) ce Peat ®) a Pe nel 


in reeller Rechnung kénnten Verallgemeinerungen der Formeln vom Abschnitt 2 
herangezogen werden. Da jedoch bei reellen Polynomen die Nullstellen von 
Syr41,...,r+m-1(*) nach der Kreisabbildung ohnehin paarweise anfallen, diirfte 
zweckmiassiger die wiederholte Abtrennung der zu je zwei festgestellten Null- 
stellen gehérigen Terme nach Abschnitt 2 sein. 

,0(x) konvergiert auch, wenn unter den m Nullstellen mehrfache sind’). 
S(z) enthalt dann ebenfalls die entsprechenden Mehrfachheiten. Auf die Bildung ~ 
von "+™Q(«) gehen wir nicht weiter ein, da dieser Fall in der numerischen Rech- — 
nung ohnehin besondere Massnahmen erfordert (vgl. § 5. 4). 


§ 4. Die abgekiirzte Iteration einer Matrix an einem Vektor 
1. Abkiirzung der gewohnlichen Vektoriteration - 


YW sei eine beliebige (7 1)-reihige Matrix mit den Eigenwerten €,(2 = 1,..., 7), © 
die so numeriert sein mégen, dass sie der Bedingung (2) gentigen. Die wieder- 
holte Anwendung der Matrix % auf einen Ausgangszeilenvektor t liefert 


+o ee Hee FoR Ce ee ae (40) 


Falls |A,| > |A,| ist, geht der normierte Vektor t,/|t;| gegen den nor-— 
mierten Linkseigenvektor u,, vorausgesetzt, dass t u, enthalt, also nicht un- 
gliicklich gewahlt ist. Dies bildet die Grundlage der Methode der Iteration an — 
einem Vektor, die von v. Mises und PoLLAczEK-GEIRINGER [11] zur numeri- 
schen Bestimmung des gréssten Eigenwertes eingefiihrt worden ist. 

Die Iteration an einem Vektor kann nun ebenso abgekiirzt werden wie die 
Iteration mod P(x). Bezeichnen wir mit T die Matrix 


to brace it, 
ho || / tu that 
T= | |» so gilt TW — | = ; . (42 42) 
ers Heats ; 
tna | tr—1 UW tign—1 
Es ist dann 


tee = t WU? = t UW T-15 = 4,5-15 ye, 


*) Vel. [4], § 4. 
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also t; 
| ties 
Tag =t tai} (43) 
ieee 
tieder || 
Kennt man also eine Folge von  Iterationsvektoren t;, ..., t;,.,_, und T-1, 


so erhalt man t,; durch gewichtete Summierung dieser Vektoren; die Gewichts- 
faktoren liefert t; I-14. t; T-1 = 3 kann natiirlich auch jeweils direkt als Lésung 
von 37 =t; erhalten werden. 

Notwendig und hinreichend fiir die Existenz von I~! ist im Falle linearer 
Elementarteiler, dass t 1 Eigenvektoren aufspannt. In jedem Fall ist notwendig, 
dass WU nicht-derogatorisch ist. 


2. Abkiirzung der Iteration an einem « abgednderten» Ausgangsvektor 


Insbesondere existiert also fiir einen neuen Ausgangsvektor "t, der aus t 
durch Annullierung der Anteile aller Eigen- und Hauptvektoren zu €,,...,&,_4 
entsteht, das zugehérige ("T)~1 nicht. Die abgekiirzte Iteration kann trotzdem 
durchgefiihrt werden, und zwar unter Verwendung der urspriinglichen Matrix 
I-1. Wir nehmen, obwohl die Beschrankung nicht notwendig ist, an, dass 


nur lineare Elementarteiler hat. Zu é,,..., €&, gehéren die Rechtseigenvektoren 
X,,-.-., ¥, und die Linkseigenvektoren 14, ..., u; eS sei (u, X,) = 1 normiert. 
Dann ist 


ae vy-1 

n= t= S'tx)y,-t(1- 5K), (44) 
1 1 

-wo K, = x, u, die Frobeniusschen Kovarianten sind, fiir die gilt 


Ko = Ke K,, K,=0 fir w+v (45) 
und 
Kt = ok, (46) 


Es sei entsprechend “t; = "t &*. Dann ist 


i OP = "t, z-1 


Tews 


Denn es ist "t;,, =+t(1— 2 K,) Wt* =t W(1 — » K,) & wegen (46). Damit 
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entsteht auf der rechten Seite 
tH Tt tad A) a 2 IC ea J) Oe Se ea 2 Ke ae 


Wegen (1— 2 K,)2=1— 2K, nach (45) ist dies aber gleich "t,,;. Falls 
\é,| >|&,41| ist, geht wiederum der normierte Vektor ’t; ;/|"t:| gegen den nor- 
mierten Eigenvektor u,. 

Es wird kaum nétig sein, zu bemerken, dass die hier beschriebene Iteration, 
mit der Transponierten %7 durchgefiihrt, auf dasselbe hinauslauft wie die 
Iteration mit 9{ an einem Spaltenvektor. Der Start eines neuen Laufes hat 
fir |é,| >|&,4| zu geschehen mit 


r+lt = 7 — (b x,) up, (48) 


erfordert also die Kenntnis des Rechts- und des Linkseigenvektors zu &, und 
damit im allgemeinen die vorherige Durchfiihrung sowohl einer Iteration an 
einem Zeilenvektor als einer solchen an einem Spaltenvektor. Fiir symmetrische 
oder antisymmetrische Matrizen WU entfallt diese Notwendigkeit. 


3. Der Fall von Eigenwerten gleichen Betrags 


Fir den Fall.) é,| = |E,:04| = += 1] éeqey | > eae) chat Are [seme 
erschépfende Diskussion der Méglichkeiten zur Bestimmung des Eigenwertes 
und der zugehérigen Eigenvektoren gegeben, und zwar auch fiir den Fall zu- 
sammenfallender Nullstellen und hdherer Elementarteiler. Wir kénnen uns 
daher mit diesem Hinweis begniigen. Fiir den Regelfall einfacher konjugiert- 
komplexer Paare von Eigenwerten und Eigenvektoren ist nach sinngemdsser 
Ubertragung der Gedankengange von § 3.2 eine rein reelle Durchfithrung der 
Rechnung méglich. 


4. Spezialfall: Die Nullstellenbestimmung von Polynomen 


Die klassische Bernoullische Methode kann bekanntlich ([1], [7]) als Vektor- 
iteration mit der Frobeniusschen begleitenden Matrix ¥ 


Diet eel 0 0 
(pee 1 ae 
ice Sas het icant ee ee (49) 
Olan 0 1 
ty =Gae4 —Qiep ess oy 


an einem Spaltenvektor (oder als Iteration mit ¥? an einem Zeilenvektor) auf- 
gefasst werden. Weniger geldufig scheint der Gedanke zu sein®), die Iteration 
mit X an einem Zeilenvektor durchzufiihren. Man erkennt unmittelbar, dass 


) Vel. [4], §3.1, Fussnote. 
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dies der in § 2 beschriebenen Iteration + mod P(x) aquivalent ist, wenn man 
den Koeffizienten des Iterationspolynoms den Iterationszeilenvektor, von 
riickwarts gelesen, zuordnet. 

Da man die transponierende Matrix von X, TX = X? T von vornherein 
kennt, 


l| @n—a Ana ay ii 
pee an—3 il 0|| 
T =| pee et aeRO \% (50) 
hie. ued ee Oe 
if 0 diet 0. 


sind beide Methoden mathematisch gleichwertig, die Zeileniteration kann in 
die Spalteniteration tibergefiihrt werden und umgekehrt. Bei der abgekiirzten 
Iteration jedoch hat die letztere, auf die Lagrangesche Variante hinauslaufende 
Methode den Vorzug, dass T = T-! = 1 wird, wenn als Anfangsvektor (dem 
Polynom Q(x) = 1 entsprechend) der Vektor (1 0 0 0... 0) gewahlt wird. Die 
aus ” sukzessiven Iterationsvektoren 


Il t, 


| ty4a 


gebildete Matrix stimmt dabei nach (42) gerade mit der 7-ten Potenz von X 
iiberein. Es ist aber bedeutsam, dass man die zweite bis 1-te Zeile von X* durch 
Iteration aus der ersten Zeile gewinnt, sowohl vom Standpunkt des Aufwands, 
verglichen mit der expliziten Potenzierung von ¥ nach AITKEN [2], als auch 
im Hinblick auf Rundungsfehler (vgl. § 5. 2). 

Es ist im tibrigen leicht nachzuweisen, dass im ersteren Fall, wenn mit ¥7 
an t = (00---01)?%) iteriert wird, T-! = T [Gleichung (50)] wird. 


§ 5. Gesichtspunkte zur praktischen Durchfiihrung des Verfahrens 


1. Gang des Rechenverfahrens 


Die abgekiirzte Iteration zur Gewinnung von Q(x) = Di ty) x-” stiitzt 
sich nach (13) und (23) auf die folgenden beiden Algorithmen (m= 1,..., 7) 


(Che) arg 


os yf) — Ay Yeo». Yoav 9, (51) 
try (B G41) zd ee ty(Bi +4) ty Barn) . (52) 


10) Diesen Fall hat schon FUrsteNAvu [8] betrachtet. 
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Ein Gesamtschritt der abgekiirzten Iteration besteht aus  Schritten (51) 
und einem darauf aufbauenden Schritt (52), der die Ausgangsdaten fiir die 
Wiederholung des Prozesses liefert. Da es auf Skalenfaktoren nicht ankommt, 
darf vor Beginn der nachsten Schritte (51) der ganze Vektor ’y{2i») durch 
eine beliebige Konstante dividiert werden, zum Beispiel durch 2er- oder 10er- 
Potenzen zwecks Einhaltung einer festen Gréssenordnung oder auch durch 
ry\8i), der Normierung von Q(x) entsprechend. 

Falls « = 0 gewahlt wird, wie es der Diskussion von § 2 entspricht, besteht 
die Méglichkeit laufender Aufsummierung aller Glieder. Die Wahl irgendeines 
t=1,...,%—1 fuhrt, wie man leicht sieht, zur Beziehung Bai =e 
erhéht also die Konvergenz geringfiigig. 

Der erste Iterationslauf kann, der Wahl Q(x) = 1 entsprechend, mit 


Oa ae 
(53) 
| 1 far» =n, 


begonnen werden. Man kann Iterationsschritte einsparen, wenn. man sogleich 
mit Q™(x) = x" — P(x) beginnt, also 


=f = -a, (54) | 


setzt. Bei einheitlicher Fiithrung aller Iterationslaufe, wie man sie bei pro- 
grammgesteuerten Rechenanlagen anstrebt, ist dann jeder neue Lauf zu be- 


ginnen mit ¥ 
7 +10)(60) (x) prs 7Q'o) ( x.) ce Eo P(x) ; (54a) 


Die Bestimmung von Dy aed, ety rem—1%) aus Ey x1, areal y ram) nach (38) 
[Spezialfalle (26) und (4)] durch K oeffizientenvergleich bedarf keiner Erlauterung. 
Mathematisch aquivalent, numerisch jedoch gelegentlich vorteilhafter ist 
es, die Koeffizienten o, von S,,, PROM (4 he he ts er ee 
aus der Relation4) " acncatliies. i a 


QM EH) +o, + 4, (QS 0 (55) 


za gewinnen. Dies kann durch Auflésung, nach Streichung iiberfliissiger 
Elemente, geschehen. Die genauesten Ergebnisse, wenn auch mit hdherem Auf- 
wand, erhalt man durch Ausgleichung, also durch Lésung des zugehérigen 
Systems von Normalgleichungen. | ; 
Bei Polynomen mit numerisch gegebenen Koeffizienten wird das Auftreten 
von Nullstellen gleichen Betrags (ausgenommen einfach konjugiert-komplexe 


sh he a ie ate ; 
) Spezialfalle (17), (27). Fir die allgemeine Vektoriteration (vgl. § 4) ist dieser Sachverhalt 


gelaufig [6], er lasst sich jedock ; : : 
elec : jedoch auch direkt aus [4], § 3, insbesondere Gleichung (66a) und (55) 
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Paare bei reellen Polynomen) meist zufallig sein. Es wird in solchen Fallen also 
geniigen, zunachst das Augenmerk nur auf die’ Konvergenz von "0 (x) oder 
allenfalls 3Q(x) zu richten und erst dann weitere, abgeleitete Folgen heranzu-’ 
ziehen, wenn jene nach einer im allgemeinen ausreichenden Zahl von Schritten 


_ nicht konvergieren. Dabei kann die Feststellung der Konvergenz, wenngleich 


weniger zuverlassig, auch am Vergleich sukzessiver Naherungswerte fiir £, oder 
P=, + 641 9=§, 6&4, erfolgen. Am genauesten ist natiirlich wiederum die 
Betrachtung des mittleren quadratischen Fehlers im zugehorigen Ausgleichs- 
problem. 

Dass auch die abgekiirzte Iteration bis auf Skalenanderungen divisionsfrei 
gefiihrt werden kann, macht sie besonders auch zur Bestimmung dei Null- 
stellen von Polynomen mit komplexen Koeffizienten geeignet. Selbstverstand- 
lich wird man bei hohen Polynomen auch daran denken, die reziproken der 
kleineren Nullstellen aus dem Polynom P-}(x) = x” P(1/x) zu bestimmen. 


2. Selbstkorrektur und stabilisierende Eigenschaft der abgekiirzten Iteration 


Die gewohnliche Iteration mit x mod P(x) hat als Spezialfall einer Vektor- 
iteration die Eigenschaft, selbstkorrigierend zu sein. Ist |é,| >|&|, so fihrt 
die Iteration fiir jedes Ausgangspolynom, das nicht &, als Nullstelle enthalt, 
zu P,(x) und €,. Obwohl unter dem Einfluss von unvermeidlichen Rundungs- 


_ fehlern die theoretische Iterationsfolge fortwahrend abgeandert wird, wiirde 


jedes neuentstehende Iterationspolynom theoretisch zur Konvergenz fihren. 
Mit Q(x) = 1 [oder auch Q(x) = P’(x)] steht ein Ausgangspolynom zur Verfii- 
gung, das im allgemeinen durch Rundungsfehler nicht mehr so weit beeinflusst 
werden kann, dass die obige Bedingung fiir die Konvergenz nach &, verletzt 
wird. Im iibrigen besteht hier sogar absolute Selbstkorrektur insofern, als es 
wegen des Zufallscharakters der Rundungsfehler unméglich ist, dass die obige 


_ Bedingung, wenn sie einmal verletzt ist, dauernd verletzt bleibt. 


Es ist verstandlich, dass die Selbstkorrektur nur méglich ist, wenn die 
«urspriingliche Information», wie es hier mit den Koeffizienten des vorgege- 
benen Polynoms der Fall ist, immer wieder in den Prozess eingreift. 

Aus der Tatsache der absoluten Selbstkorrektur geht deutlich hervor, dass 
die Iteration mit x mod P(x), begonnen an einem abgednderten ’Q(x), die 
theoretisch fiir |&,| >|é,,,| zu B(x) und &, fiihren sollte, unter dem Einfluss 


der Rundungsfehler in die Iterationsfolge nach P,(x) und &, einmiindet, also 


«labil» ist. Als stabilisierend bezeichnen wir eine zusdtzliche Massnahme, durch 


die ein Verfahren auf einen Weg, den es unter dem Einfluss von Rundungs- 


fehlern verlassen méchte, gezwungen wird. Solche Massnahmen, wie sie fiir 


Vektoriteration etwa das Verfahren von Kocu beniitzt, sind meist lastig. Es ist 
nun ein besonders hiibscher Zug der Abkiirzung, die die quadratische Konver- 


-genz liefert, dass sie gleichzeitig stabilisierend wirkt. Selbstverstandlich kann 


ek ak ee 
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das nur dadurch geschehen, dass die Verscharfung zur absoluten Selbstkorrek- 
tur preisgegeben wird. In der Praxis ist das bedeutungslos, da mit oder ohne 
diese Eigenschaft grosse Verfalschungen, wie sie durch Rechenirrtiimer entstehen, 
so spiirbare Riickschlage bedeuten, dass sie durch Verwendung von Rechen- 
proben verhindert werden miissen. Wir werden aber zeigen, dass die abgektirzte 
Iteration selbstkorrigierend bleibt, so dass der Prozess fiir ein Ausgangspolynom 
'O(x) auf P(x) und &, zielt. Wir untersuchen dazu die Fortpflanzung einer 
, Verfalschung bei der abgekiirzten Iteration. Das Iterationspolynom 


10 (x) = ar E Pls) 


sei durch Rundungsfehler entstellt, statt dessen liegt das Polynom 


00) = Ss, Pla) + S04 o) BBW) (6) 


vor, wo |e, &!| <1 ist. Wird nun durch Iteration mit x mod P(x) die theore- 
tische Folge ’Q*(x),...,’Q4*"-(x) hergestellt, so ergibt sich durch ge- 
wichtete Summierung 


y—1 


Orns) = Se)? P(x) + V+)? B(w). (67) 


Man erkennt die stabilisierende Eigenschaft darin, dass die Anzahl von Stellen, 
die von Stérkomponenten in den ersten v Teilrdumen nicht beeinflusst werden, 
nach jedem Gesamtschritt etwa doppelt so gross ist als vorher und praktisch 
der Stéreinfluss (bis auf neu entstehende Rundungsfehler) jeweils eliminiert wird. 
Nur bei sehr nahe benachbarten Nullstellen und dann am Anfang der 
Iteration kann das Verfahren in eine andere Lésung ausbrechen [wenn etwa 
(1+ ¢,) &< (1+ ¢,,,) &,, wird!. Aber selbst dies bedeutet nur einen Wechsel 
in der Reihenfolge, in der die Nullstellen bestimmt werden. Ein labiler Zustand, 
etwa das Gleichheitszeichen in der vorstehenden Ungleichung, kann unter dem 
Einfluss von Rundungsfehlern wiederum auf die Dauer nicht bestehen. 


3. Rechenproben 


Fir die Iteration mit x mod P(x) existiert eine Rechenprobe, die mit einer 
Summenprobe der entsprechenden Matrixiteration zusammenfallt, naémlich 


Qh) = 08) = ye Pa): (58) 


ebenso hat man fiir die gewichtete Summierung die Summenprobe 


7921) = DS: Foy) 16 a= (iN (59) 
1 
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4. Stellenverlust durch Ausléschung 


Ein Verlust an genauen Stellen durch Ausloschung tritt beim abkiirzenden 
Schritt (52) nur ein, wenn der zu bestimmenden Nullstelle &, eine oder mehrere 
Nullstellen relativ zu ihrem Betrag nahe benachbart sind. Man erkennt dies aus 
der Betrachtung des (fiir numerisch gegebene Polynome meist idealisierten) 
Grenzfalls mehrfacher Nullstellen?”). Es treten dann hohere Elementarteiler der 
begleitenden Matrix X auf; konvergenzbeschleunigende Prozesse, wie die Ab- 
ktirzung oder die Bildung der ersten normal konvergenten Folge unter den ab- 
geleiteten Iterationsfolgen, fiihren zwangslaufig zum Auftreten von Differenzen 
annahernd gleicher Zahlen, also zu einem Genauigkeitsverlust durch Ausléschung 
fihrender Ziffern. 

Bei der Wahl von P’(x) als Ausgangspolynom wird das Auftreten von 
Elementarteilern theoretisch zwar verhindert, praktisch jedoch im Hinblick 
auf Rundungsfehler nur ein labiler Zustand hergestellt. Die Abkiirzung der 
Iteration als Spezialfall von § 4 ist in diesem Fall — wenigstens unmittelbar — 
nicht méglich, da die zu verwendende Matrix T~—1! gerade im Falle mehrfacher 
Nullstellen nicht existiert. Es besteht jedoch eine Méglichkeit, die iiber den 
bisher diskutierten Rahmen hinausgeht. Wird namlich der unveranderte Pro- 
zess der abgekiirzten Iteration von § 2, beginnend mit dem Ausgangspolynom 
Q':)(x) = P’(x) oder allgemeiner = x* P’(x) mod P(x), durchgefiihrt, so ent- 
scheidet tiber die Konvergenz, welcher der Ausdriicke &, P’(&,,) (m= 1,..., ”) 
betragsmassig am gréssten ist. Da P’(é,,) fiir relativ nah benachbarte Null- 
stellen sehr klein ist und fiir mehrfache ganz verschwindet, greift das Verfahren 
nun zunachst die am besten isolierten Nullstellen an und zeigt dabei tiberhaupt 
keine Elementarteilerschwierigkeiten. 


5. Eigenwertproblem 


Grundsiatzlich iibertragt sich die Diskussion der Abschnitte 1 bis 4 sinn- 
gemass auf die praktische Lésung von Eigenwertproblemen nach § 4, insbe- 
sondere hat man nach (47) wiederum eine Stabilisierung durch den Abkiir- 
zungsschritt. Wenn neben den Eigenwerten auch die zugehérigen Rechts- und 
Linkseigenvektoren gesucht sind, hat das Verfahren den Vorteil, dass wegen 
der Selbstkorrektur eine Nachiteration entfallen kann. Aufwandsmassig ist 
jedoch fast immer der Weg iiber die Bestimmung des Sakularpolynoms giin- 


stiger, da man dieses schon nach dem Verfahren von FRAZER, DUNCAN und 
n-1 
CoLaR [9] als Lésung von »’ a, t, = —t, mit demselben Aufwand erhalt, den 
d 0 


| Bac einzige Berechnung von ts, T~' erfordert. Zur Bestimmung der Eigenvek- 
toren ist die unmittelbare Kenntnis der Interpolationspolynome P,(x) von 


#2) Vel. [4], § 4. 
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besonderem Vorteil; im Falle einfacher Eigenwerte hat man bekanntlich 


Po) =a ates: (60) 


(6 


Man kann diese Bezichung auf einen Zeilen- oder Spaltenvektor anwenden, 
um u, bzw. x, zu erhalten. Wiederholt man den Prozess, so erhalt man damit 


auch eine wirksame Nachiteration. 
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Summary 


One complete step of the abbreviated iteration consists of 1 single steps 
(n being the degree of the polynomial) of the common Bernoullian iteration in a 
variant based on LAGRANGE’S method. This kind of iteration, the advantages of 
which seem to be little known, particularly makes possible the aforementioned 
abbreviation in the form of a weighted summation, the weight factors being 
determined immediately during the process. 

But for trivial changes of scale factors the entire iteration may be set up 
without divisions. Therefore it is well adapted to the determination of zeros of 
polynomials with complex coefficients. 

Our iteration process gives not only the zero to be determined but also directly 
the residual polynomial without the respective linear factor. In the case of two or 
more dominating roots of equal modulus the iteration process does not converge, 
but by successive processes analogous to the ‘cross-multiplication’ of RouTH a 
derived sequence may be formed which again converges to the polynomial of the 
remaining zeros. ; 

The original polynomial may be used without any change for the determina- 


tion of each zero and there is no carrying of rounding error. Each iteration run is 
self-correcting and may be checked easily. 


(Eingegangen: 6. Dezember 1954.) 
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On the Error Function-of a Complex Argument 


By JoserH KEstTin, Providence, R. I., USA), and 
Leir N. Persen, Trondheim, Norway?) 


Introduction 


The error function of a complex argument and allied functions occur in 
many problems of physics or engineering. They occur, for example, in the 
theory of heat conduction, particularly when non-steady heat conduction 
through stratified media is concerned as well as in problems of electro-chemical 
diffusion. Among the additional fields of application, electrical transients, 
electromagnetic theory, and rocket ballistics may be mentioned. The application 
of the error function to problems in statistics and of the allied Fresnel integrals 
to classical diffraction problems is also very well known. 

In spite of their great usefulness, the tabulation of the error and related 
functions of a complex argument has not progressed very far because of their 
complexity as well as because of the large amount of numerical work involved. 
Rosser {1]8) made a comprehensive study of the methods of integrating the 
functions 


Zz loc) 


[exp (—&?) dé and exp (2%) [ exp (#9) dé 
0 


z 


and prepared suitable auxiliary tables. CLEmmMow and Miss MUNFORD [2] gave 
a short table of values of the function 


G(z) = exp( in?) / exp (—i £2) dé 
2Va/2 


[5 mesol 5 fe) foxo(—3 fa) a 


_ SKWIRZYNSKI [3] prepared a short table of the complementary error function 


I 


2% 
Ls 


erfcz = Z_ | exp (-€) Gee” Pep cide 
Vx: V 100 


1) Brown University. 
2) Norges tekniske hégskole. 
3) Numbers in brackets refer to References, page 39. 
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In all these papers the integration has been performed by more or less direct 
methods, principally by a series expansion. Recently CAHILL of the National 
Bureau of Standards [4] has undertaken a programme of computation with the 
aid of SEAC by the general code for evaluating confluent hypergeometric 
functions, making use of the identity 


© Be il Seen ei 
erfz = a A(> > 2 ) 2 ) . 
Finally it should be mentioned that LaiBLE [5] prepared a relief chart of 


the error function P 


erfz = ja] ee) dé , 


M4 
0 


calculated the first five zeros, gave an asymptotic expression for the zeros, 
and indicated a graphical method for the tracing of lines of constant modulus 
and phase angle. 
The preceding methods, apart from being tedious and time-consuming, do 
not lead to simple expressions for the real and imaginary parts of the function. 
In the present note we wish to give a transformation which seems to us to 
lead to a simple expression for the related function 


K (2) = exp (z2)rerie(2)e (14 


Explicit formulae for the real and imaginary parts of this function 
K(z) = K,(x, y) +1 K;(x, y) (2) 
Z=x+1Yy (2a) 


are obtained in terms of elementary functions and two definite integrals. The 
resulting formulae facilitate integration along lines passing through the origin, 
i.e. for z varying at a constant phase angle. 


The Transformation 


Recently the authors encountered a problem involving the calculation of 
slow torsional oscillations of a disk [6] and of bodies of revolution [7] where 
the transformation now to be explained rendered good service. 


Let « and « denote two complex conjugate numbers and consider the two 
identities: 


ei ies ao 
Vs+a sta (Vs+a) (Ys+a) ’ | 
Bou oe AR eee (ee | o 


Vota “Vs shay Vere ty eee 
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Putting 


A=a+tb, a=ayrtd, (4) 


we can rewrite the identities (3) to yield 


Dig pee as ae pea 21 8 

Vsta  ys+e (Vs + a) (Vs +a) ’ | 
1 1 AY sp 2 6) 
S(s) = — ape gt en i Lately | 

ysta Vs+e (V/s + a) (/s + a) 


We now take the inverse Laplace transforms of the left- and right-hand 
sides of equations (5). The left-hand sides, denoted by the subscript J can be 
integrated by reference to tables [8], [9], and we have 


go rae 7 ma a FlaVt) , (6) 


where 


n 


P (z) =exp(2*) eric (z) . (7) 
Employing the notation 
L1H Dis) }= Diy) , ete: , 


we have 


I 


D,(i) =aF(aVt)—«F(aVi), St) =—-—aF(aVt)+aF(aVi). (8) 


The right-hand sides of the identities in equations (5) can now be integrated 
directly noting that the integrands have simple poles at s,= «? and s, = a 
and branch points at the origin. The residues will be expressed by exponentials, 
and hence will lead to damped harmonic oscillations. The contributions from 
the branch points will be given in the form of definite integrals, as shown in 
the Appendix. 
Performing the somewhat lengthy but standard computations, we find 


K,(a Vt) a a +2yt exp| (a? — 6*) t]{ a cos(2 a bt) — bsin(2a bt)} a 
+ (b2 — a®) Vt P,(a, b, t) — Vt P,(a 62) 


and 
K (a2) =-+ 2 t exp| (a? — b) t]{ asin (2a 6t) + bcos (2a bi) }| 


: (9a) 
+2abVtP,(a, b,t), 
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where 00 a ‘ 
P,(a, b, t) a =a ee oe 2 aE am b2)2 , (10) 
0 
oo wee bee 
P,(a, b, t) = — ‘| os aS ne =a (10a) 


all the quantities now being real. It should be noted that the integrands in 
equations (10) and (10a) are always positive. The calculation leads, of course, also 


to an expression for K(« Vt) which, however, need not be considered as it follows 
from K(« Vt) by substituting & for «. 


Alternative Forms 
In order to obtain formulae for K(z), with z= x +7, we put 
b= tet VieWjamece (11) 


so that x+1y=a Vt. Hence 


K,(%, y) = sere exp (x? — y?){ x cos(2 x y) — ysin{2 x y) } | 
ya (12) 
+ +=" Pix, 9) — 9 Paley) | 
K,(x, y) = + 2exp(x? — y?) { x sin (2 x y) + y cos(2 x y) } 
+ (12a) 
+ 2% P(x, 9), 
a a f 1? exp(—€ y?) dé 
ced ae / STF 2 (aly) — 1 E+ Ua? + UE” ea 
it ~ £312 emety (eae AD UTE 
P,(x, as CRee ee S Di=s y*) dé or 
NG y) a & Tee) [(x]/y)2 aS 1] E ae [(a/y)2 Le Ta (12c) 


It is, of course, sufficient to perform calculations in the first quadrant, as 
the representations in the remaining quadrants are mirror images with respect 
to the axes of co-ordinates. The zeros of K(z), which are identical with those 
of erfc z are given by solving the two simultaneous equations 


K,(x,y)=0, K,{x, y)=0. (13) 


For large values of ¢, i.e. for large values of r = Vx2 + 2 the terms with the 
functions P, and P, may be neglected. Hence the asymptotic zeros are found 


a oe i EEE EE oO lle 
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from the system 


x cos (2 x y) — ysin(2 x y) = — ,  tan(2 x y) = 2. (13a) 
The forms (12) to (12c) give a better insight into the behaviour of the 
function A(z) and hence into that of F(z) and erfc(z) than the standard 
formulae because the oscillation appears explicitly and not under the integral 
sign. Corresponding expressions for the real and imaginary parts of F(z), 
erfc(z), and erf(z) can be easily deduced. 
Expressions in polar co-ordinates are obtained by putting 


z=re'®, ie. a=cosy; b=sing; Vi=r. (14) 
Hence 
Kor, @) = = + 2rexp(r? cos2 ¢) cos(p + y) 
ya (14a) 


—4rP,(r, y) — r(cos2 g) P,(r, g) , 
K,(r, gp) = + 2rexp(r? cos2 g) sin(y + py) + r(sin2 g) Py(7, py), (14b) 


where 
pir sin 2g ; (14c) 
EL) ( )d 
i ARES exp (—r? &) d& 14d 
Pir, @) ==] £9 2 F ces? pt 1% ee) 
0 
if &/ exp(—r? £) dé 
P,(r, 9) = | £2 bcos? oa 1 CA? 


0 


The Fresnel integral is related to the value of K(z) for y = 2/4. In this 
particular case 


Sac: ‘ 1 / 8 exp(—r? &) dé 
K,(r, z) = aie eth cos( 7 at 7?) hs =| e241 : 
* 0 
; 1 f ét exp(—r? &) dé 
Kr, *) = 4 2rsin(S +r) +7 241 ‘ 
0 
APPENDIX 


Derivation of the Result Used in Equations (9) and (9a) 


We now proceed to sketch the method of finding the inverse Laplace 
~ transforms of the right-hand sides of equations (3), i.e. to calculate D, and S,.. 
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From the first of equations (5), we have 


1 Caves 2i1b est ds 
tes a A-1 
D,(8) B i f (/s + «) (V/s + &) Ma 


The integrand has simple poles at s, = «? and s, = « and a branch point at 
the origin. The residues are: 


at OSpe=0" 
R(s,) = —2aexp (a? ?), (A-2) 
at So =a" 


R(sy) = +20 exp(x*?) . (A-3) 


Taking a cut along the real negative axis in the s-plane, and putting 
s = &e?®, we find that : 


(Vs +a) (Vs +%) =£exp(i O) + («+ %) éI” exp(5 i) + au 
for 0 = —x 
(Vs+a) (Vs+a)=-&-(0+a) fH itaa, 
forO = +2: 
(Vs+a) (Vi tole Pee eee ae, 


and the contribution from the branch point becomes 


joule. aoe xi} 2ibexp(—&t) d& 
: ($+ aa) -—i@+ = 


[oe] 
1 [- 21bexp(—é?#) dé 
( a 


ce Ee —&+aa)+i( 


Evaluating and substituting « =a+ibanda—a—i b, we obtain 


(A+) 


Thus 


D,(t) = 2% exp (&?t) — 2% exp (a2 t) — sp 2 exp(—€#) dé 
&2 4 2 (a2 — B®) & (a? + b2)2 


(A-5) 


_ 
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_ Substituting for « and % we can transform this to 
D,(t) = — 47{ exp[(a? — b) é][asin(2 abt) +5 cos(2abt)])+abP,}, (A-6) 


where P, was given in equation (10), 
Proceeding in a similar manner, we find 


si) = 2. [ (2VS+ 2a) exp(s#) ds ; 
: raf Ys+a)Ys+a) | a 


Making use of the preceding results, it can be shown that 
S,(t) = —4 exp[(a? — 6?) ¢][a cos(2 a bt) — bsin (2 a b2)] 
A-8 
—2 (b?— a) P,+2P,, | a 


where P,(a, b, ¢) was given in equation (10a). Equating D, = D, and S, = S, 
we are led to the result in equations (9) and (9a). 
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Zusammenfassung 


In der Einleitung werden die wichtigsten bisherigen Arbeiten tiber die Fehler- 
funktion mit komplexem Argument kurz besprochen. Die der vorliegenden Unter- 
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suchung unterworfene Funktion A(z) ist in Gleichung (1) definiert. Mit Hilfe der 
Laplace-Transformation und ihrer Umkehrformel wird gezeigt, wie sich die zwei 
Identitaten in Gleichung (3) behandeln lassen, so dass sowohl der Real- als auch 
der Imaginarteil der Funktion K(z) sich in zwei Teile aufspalten lasst, was aus 
den Gleichungen (9) und (9a) hervorgeht. Der erste Teil kann durch elementare 
Funktionen ausgedriickt werden, waihrend der zweite Teil sich durch zwei 
Integrale darstellen lasst. Fiir die praktische Anwendung der angegebenen Aus- 
driicke sind in den Gleichungen (12), (12a) und (14a), (14b) die Ausdriicke in 
kartesischen bzw. polaren Koordinaten umgeschrieben worden. Der Vorteil der 
angegebenen Aufspaltung liegt darin, dass die in den Ausdriicken auftretenden 
Integrale monoton abnehmende Funktionen der unabhangigen Veranderlichen 
(x, y) darstellen und sich deswegen leicht numerisch ausrechnen lassen. Der 
Schwingungsanteil der Funktion K(z) ist ausschliesslich durch elementare Funk- 
tionen ausgedriickt. Im Appendix ist der Rechnungsvorgang, der zu den ange- 
gebenen Ausdriicken fiihrt, naher umschrieben. 


(Received: February 10, 1955.) 


Aquivalenzsatz, Ahnlichkeitssatze 
fir schallnahe Geschwindigkeiten und Widerstand 
nicht angestellter Korper kleiner Spannweite 


Von FRIEDRICH KEUNE und Kraus OswatitscH, Stockholm?) 


Die erste Arbeit iiber die Ubertragbarkeit von theoretischen oder experi- 
mentellen Resultaten an nicht angestellten Rotationskérpern auf Kérper 
kleiner Spannweite stammt von Warp [1]?). Er bewies mit Hilfe der Laplace- 
Transformationen und des Impulssatzes, dass der Widerstand eines Korpers 
kleiner Spannweite in Uberschallstrémung innerhalb des Linearisierungsberei- 
ches der Differentialgleichungen gleich ist dem:+Widerstand eines Rotations- 
kérpers gleicher Querschnittverteilung Q(x). Solche Kérper gleicher Verteilung 
der Querschnittflachen-Inhalte (Fig. 1) sind der Angelpunkt aller folgenden 
Uberlegungen. Sie seien als «Aquivalente Kérper» kleiner Spannweite’) be- 
zeichnet. Eine wesentliche Voraussetzung des Satzes von WARD liegt in der 
Forderung, dass die Querschnittanderung am Kérperende x = X» verschwinden 
muss [Q,(%9) = 0], das heisst der Aquivalente Rotationskérper muss entweder 
zylindrisch oder kegelig enden (Figur 2a und 2c). | 

Unabhangig von WarD und auf anderem Wege wurden Korper kleiner 
Spannweite im Linearisierungsbereich der Unter- und Uberschallstroémung von 


1) Kungl. Tekniska Hégskolan, Institution Gasdynamik. 


°) Die Ziffern in eckigen Klammern [ ] verwei f i ichni i 
ieee haduet aces il isen auf das Literaturverzeichnis auf Seite 63. 
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KEUNE und OswatitTscu [2], [3] untersucht und die Rechnungen auf die Ge- 
schwindigkeitsverteilung solcher Kérper ausgedehnt. Ausgehend von den Quell- 
belegungsintegralen fiir diinne Fliigel kleiner Spannweite*) wurde gezeigt, dass 
sich die Stérung der w-Komponente auf zwei Effekte zuriickfiihren lasst. Der 


a) 


Rb) 


Ox, )> oO 
Af) =15} 


b) 


Wt CPx) — x 
Q(x) >O 
Qy [Xq)#O 


2 Rl) 


xX 


Wig) = lr)=0 


at, @ 
%*o 

Qlxo)=0 

Qf x, ) Ee o 


Figur 1 


Fliigel kleiner Spannweite und Aquivalenter Figur 2 
Rotationsk6rper. Die vier Typen von Rotationskérpern. 


eine Anteil, die «zweidimensionale Querschnittstrémung», gehorcht der La- 
place-Gleichung fiir das Stérungspotential (co Index fiir das Anstrémgebiet) 


Go 1=0% Go =% Goth (1) 


Ue Ugo Uo 


in einzelnen Querschnitten 
Puy + Pez = O- (2) 


Dieser Teil erfiillt die Randbedingung am Korper in den einzelnen Quer- 
schnittebenen. Er ist fiir Uberschall- und Unterschallstrémung der gleiche. 


4) Siehe etwa [4], S. 373 ff. und S. 383 ff. 
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Beim Problem angestellter unendlich diinner Platten, also in der hier nicht 
behandelten einfachen Theorie von JoNEs [5] tritt iberhaupt nur diese Quer- 
schnittstrémung auf. Sie fiihrt auch in unmittelbarer Schallnahe zum Resultat, 
denn in Kérpernihe lasst sich die linearisierte wie die nichtlinearisierte gas- 
dynamische Gleichung in diesem Fall auf die einfache Form (2) reduzieren. 

Zu dieser die Kérperform beriicksichtigenden « Querschnittstromung» tritt 
bei nichtangestellten Kérpern noch ein zweiter Effekt, der « Raumeinfluss». Er 
ergibt sich aus den Anderungen des Kérperquerschnittes in Stromungsrichtung. 
Da seine Beitrige von den entlegeneren Kérperteilen herstammen, hat er achsen- 
symmetrischen Charakter. Es zeigt sich, dass dieser Raumeinfluss fiir alle 
Aquivalenten Korper gleich ist. 


Verkeiling der Querschnitt -Flacheninhalle 


Rotationskorper: os egal e (8) 
(A)-):Grndrissy =O al 
/5) | 
2¥i | f ) 88: | 
Hi | ce eo . 
Pas ) a | 
3 | 
+10 \ zp 
Hy mie 
ifaw (8)-1B)Mar Querschnit Fy 0) 
! xnx, 4 (c) 
i 6 ee 
(A) == Z) 
Ha p> ; : = ha 
LG 
(ce) . (€)-(c): Langsscinilt:z =O (c) 
| (8) a) 
+40. | 
Rix) (A) 
6 (8) He if (B ie 
10 jl dal ) 
4) ea) 
-10 %s ae: (8) ae 
8) pe Sg o () 
10 4 ae! 
(8) 
Figur 3 


Reispiel BE ccatala re merper. Der Spindel als Aquivalentem Rotationskérper [O(@) =146 

x *« (1—* ) mit der K6rpertiefe als Einheit] sind drei Fligel a rie G ad i ee 

eee aes Die Querschnittform ist zunachst beliebig. Die Seat des Fliigels in 

cs, glen aeons von der Form eines Parabelbogen-Zweieckes gezeichnet. Q pe b 
- Pur 6 = 0,40 und Q,,, = 0,02 sind die Langsschnitte unverzerrt Bae. 


EE EE : 
EEE eee 
—— 


ee ee 
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Am 8. Internationalen Kongress fiir theoretische und angewandte Mechanik 
hat OswatitscH gezeigt [6], dass innerhalb einer angebbaren Genauigkeits- 
grenze der Satz von der Gleichheit der Raumeinfliisse aquivalenter K6rper 
auch in Schallnahe gilt. Dieser Aquivalenzsatz soll zunichst in anschaulicher 
Weise hergeleitet werden. Sein exakter, an mehreren Beispielen erlduterter 
Beweis erschien an anderer Stelle [7]. Figur 3 gibt einen Einblick in die Mannig-- 
faltigkeit geometrischer Formen «dquivalenter Kérper». 

Mit Hilfe des Aquivalenzsatzes ist die Geschwindigkeitsverteilung an Delta- 
fliigeln®), Pfeilfliigeln und ahnlichem mit den in Schallnahe tiblichen Spann- 
weiten auf die Geschwindigkeitsverteilung an Rotationskérpern zuriickgefiihrt. 
Die Bestimmung der letzteren bleibt einem allerdings nicht erspart. Mit Hilfe 
von Ahnlichkeitssatzen fiir Rotationskérper [8] und des Aquivalenzsatzes ist 
es jedoch méglich, theoretische oder experimentelle Ergebnisse auf Fliigel 
oder K6érper umzurechnen, welche durch affine Verzerrung in Spannweiten- 
oder Dickenrichtung auseinander hervorgehen. 


1. Der Aquivalenzsatz 


Der K6rper (Figur 1) sei im folgenden durch Quellbelegung in seiner 
Grundrissflache [(x, z)-Ebene, z in Spannweitenrichtung] gegeben: 


atx 2 oe, (3) 


ten 


dabei sind die Komponenten wu, v, w der Stérgeschwindigkeit parallel zu den 
Achsen x, y, z angenommen. Mit der halben lokalen Fliigeldicke 


y = h(x, 2) (4) 


ist die 6rtliche Quellstarke bei diinnen Fliigeln direkt durch die Dickenanderung 
in Anstrémrichtung (x) bestimmt: 


q(x, 2) = 2h,(x, 2) - (5) 


Fiir die Anderung der Querschnittflache Q,,(x) gilt fiir diinne und dicke Fliigel 
+a(x) 
Qalx) = [ a(t) dt. (6) 
—oa(*x) 


Die halbe Breite o(x) der Quellbelegung in Spannweitenrichtung stimmt bei 
diinnen Fliigeln mit der lokalen Halbspannweite tiberein. Fir dicke Fligel mag 


5) Unter « Deltafliigel» werden hier nicht nur Dreieckfligel, sondern alle die Fliigel verstanden, 
die in der Vorderkante keine endliche Spannweite haben («pointed bodies»). 


44 FRIEDRICH KEUNE und KLAus OSwATITSCH ZAMP 


sie grésser sein als o(x), ausserdem ist dann Gleichung (5) nicht mehr genau 
genug®). Wie Kérperform und g(x, #) in diesen Fallen zusammenhangen, braucht 
hier jedoch nicht untersucht zu werden. Fir jeden Korper kleiner Spannweite 
ist der Aquivalente Rotationsk6érper durch die Quellbelegung mit Gleichung (6) 
gegeben. 

Aus der Quellbelegung folgt analog zur Profiltheorie fiir das zweidimensio- 
nale Potential der Querschnittstroémung an jeder Stelle ~ des Korpers: 


+a(%) i Ly WD 
p(x, y; 2) Ti = f q(x, t) In) y2 AR (t i 2)? dt 5 (7) 
— a(x) 


Es erfiillt die Laplace-Gleichung (2) und liefert die nach Gleichung (3) ge- 
forderten v-Komponenten auf der Projektionsebene. Die u-Komponente der 
Querschnittstr6mung ist dann 

: a 5(2) “ 
—oa(x) : 


Die Differentiation ist wegen der von x abhangigen Grenzen nicht unter das 
Integral gezogen. 

Beim aquivalenten Rotationskérper konzentriert sich die Quellbelegung 
auf einen schmalen Streifen 2 ¢ an der Achse, und man erhalt mit Gleichung (6): 


+e 
1 TASER pe OO, i erie ees 
Pree Kee a In y2 tae Ae [ae t) Gta A In)'y? cig Se (9) 


2a 


Der Raumeinfluss der u-Komponente ist definitionsgemass die Differenz 
zwischen gesamter u-Stérung und ebener Querschnittstromung. Da der Raum- 
einfluss in Kérpernahe unabhangig vom Achsenabstand ist, kann sein Wert 


auf der Fliigeloberflache y= h mit dem Wert auf der Oberfliche des Aqui- 
valenten Rotationskérpers 


ee ee eee 


Vy?+22= R(x), O(x) <2 R? (10) : 


verglichen werden. Im folgenden sollen die Werte auf der Projektionsflache des 
Fliigels (y = 0) betrachtet werden, die bei diinnen Fligeln mit jenen auf der 
Oberflache nahezu iibereinstimmen. Bei dicken Fliigeln sind die Unterschiede 
zwischen den Werten auf der Grundrissflache und auf der Oberflache des 
Fliigels nach KEuNE [9], [10] durch einfache Beziehungen bestimmt. 


6) Vergleiche [10], Gleichung 14. 
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Die Definition des Raumeinflusses fiihrt damit unter Verwendung von 


Gleichung (8) und Gleichung (9) zur analytischen Formulierung des Aquivalenz- 
satzes’) : 


a(x) 
u(x, 0, 2) 1 a> A 
( = 1) mae rs ‘ h,(x, t) In|t — 2| dé 


— a(x) 


A) oe ew Goel) RG). (A) 
Zu einer Fehlerabschatzung fiir Gleichung (11) fiihrt die folgende Uberle- 
gung. 
Die ¢,-Verteilung des Fliigels, Gleichung (8), geht bereits in wenigen o(x)- 
Abstanden von der x-Achse sehr rasch in die achsensymmetrische Verteilung, 
Gleichung (9), iiber, wie eine Entwicklung des Logarithmus fiir 


Vy +22>o0(x) St 
zeigt. 

Daraus folgt aber die einleuchtende Tatsache, dass sich die Geschwindig- 
keitsverteilung des Fliigels kleiner Spannweite von jener des Rotationskérpers 
nur innerhalb eines koaxialen Zylinders wesentlich unterscheidet, der einen 
Durchmesser von der Gréssenordnung der Spannweite 2) hat (Figur 4). 


Figur 4 


Bereich der Kompressibilitatsquellen. 


Wird die vereinfachte gasdynamische Gleichung 


TOA ie rd ree eg oma 


Ox Uno Uoo Uoo 


fiir den Raumeinfluss geschrieben, so erhalt sie mit Riicksicht auf Gleichung (2) 


7) Vergleiche [7], Gleichung 75. 
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die Form: : 
a—m) x [(* -1)- 9] + oy Lee ~ MI : 
+4; ies p.| = —(1— M) ore 3) 


mit einem « Quellglied» auf der rechten Seite. 

Gleichung (13) gilt sowohl fiir den Kérper kleiner Spannweite als auch fiir 
den aquivalenten Rotationskérper, indem auf der rechten Seite von Gleichung 
(13) @z2 nach Gleichung (8) und Gleichung (9) zu setzen ist. . 

Im Fall der linearen Gleichung mit 1— M*= const fiihrt die Differenz 


beider Gleichungen zur gasdynamischen Gleichung fiir die Differenz der ee 
einfliisse mit einem Quellglied auf der rechten Seite, welches nur in dem in 
Figur 4 abgegrenzten Raume wesentlich in Erscheinung tritt. Dabei sind die! 
Beitrage in Achsennahe wegen ihrer geringen raumlichen Ausdehnung beson- 
ders gering. Bei einer Berechnung sind die y- und die z-Achse am besten wie 
bei der Prandtl-Glauert-Analogie affin zu verzerren. Daraus folgt dann als 
Quellbeitrag die Gréssenordnung von gy, des Fliigels kleiner Spannweite 


Gael) Wy)5/1 — Me) 9), 
multipliziert mit dem um den Faktor 1 — M? vergrésserten Volumen (vgl. 
Figur 4 mit x) = 1). So kommt man im linearen Fall zum Fehler 

b2 (1 — M2) eae In|b/1— M2| (14) 
von Gleichung (11). 

Im nichtlinearen, schallnahen Gebiet bleibt die Tatsache bestehen, dass die 
wesentlichsten Unterschiede in den beiden zu vergleichenden Strémungsformen 
in der Nahe der Oberflache des 4quivalenten Rotationskérpers auftreten. Wird 
die Strémung am Fliigel kleiner Spannweite als Stérung der Str6mung am 
Rotationskérper aufgefasst, so leuchtet es ein, dass bei der Entwicklung der 
Unterschiede beider Strémungen die mittlere Machzahl in K6érpernahe des 
aquivalenten Rotationskérpers ausschlaggebend sein muss. Das heisst jedoch, 
dass mit einer solchen Machzahl die Abschatzung Gleichung (14) auch fiir die 
gesamte Schallniahe gilt. 


2. Schallnahe Ahnlichkeit 


; Der Ableitung der Ahnlichkeitsgesetze fiir Fliigel klemer Spannweite sei 
eine kurze Ableitung der schallnahen Ahnlichkeit fiir Rotationskérper [8] voran- 


8) Mit Q,_(0)/2% als Quadrat des Tangens des halben 


Offn: i Aqui 
Rotationsk6rpers, Gleichung (19). ungswinkels des dquivalente 
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gestellt. Wie in zahlreichen fritheren Arbeiten®) sei der Faktor 1 — M2 durch 
eine lineare Funktion von w genahert: 
aa 2 7 pa 2 = faa Uoo 
1a = (1 2M) [2 a a (15) 

c* ist die kritische Geschwindigkeit, M,, und u,, der Anstrémzustand. Man 
beachte, dass Gleichung (15) auch fiir M,, > 1 sinnvoll bleibt, weil gleichzeitig 
t.. > c* gilt. 

Die gasdynamische Gleichung (12) und die Gleichung fiir Wirbelfreiheit 
(y? = y? + 2?) lauten dann mit Gleichung (15): 


r[a M?,) poke (“= *e) | 0 (<*) =0; 


Tico, Ox Uo 


16 
—s Ub — Ups 0 af (16) 
wo aS ek “ue : 
wenn 
eo, SSS pero : 
A (ij) —1 wnetl (17) 


ist. Im Unterschall miissen in grossem Abstand vom Korper alle Stérungen 
abgeklungen sein, eine hier stets erfiillte Randbedingung. Die Randbedingung 
auf der Achse ist: 

Ae} Gakgae TF ne (18) 


Uoo Pasi 


Zur Einfiihrung neuer dimensionsloser Grdéssen wird die geometrische 
Grésse Q,,,(0)/2 herangezogen, welche mit dem halben Offnungswinkel 
der Kérperspitze des Aquivalenten Rotationskérpers wie folgt zusammenhangt: 


Beal’). — tg, (19) 

Dann seien: 
rea aire nies Hos, RE oo el ae ne Ba 
‘dara a? u ie 0,,(0) vr vie 0,,,(0) ( ) 


Die neuen Ausdriicke fiir die Geschwindigkeiten liegen im Wertebereich der 
Einheit. In den Variablen (20) lauten die Gleichungen (16) und (18): 


.[ 2a Sa EL MN 
‘Loy Me) al apt-get 
; ~ OU OD ae 
7 i, Ayer 
sidigilecm ti Qalt) (22) 
feals wii bane Rael y 


9) Vergleiche [4], Gleichung (IX, 9). 
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Affin verdickte Kérper haben alle die gleiche Randbedingung (22). Ahnliche 
Lésungen erhalt man aus der zusitzlichen Forderung der Gleichheit der Diffe- 
rentialgleichungen (21). 

Eine Verschiedenheit der spezifischen Warmen x und der Werte A? in den 
Vergleichsbeispielen kénnte durch Ubergang zu den Gréssen u A? und v A? in 
die Ahnlichkeitsgesetze einbezogen werden. Da diese weniger wichtige Verall- 
gemeinerung die Ableitung belastet, wird auf sie verzichtet, und sowohl x als 
auch A? werden als gleichbleibend vorausgesetzt. 

Die Ubereinstimmung der Differentialgleichungen (21) ist gewdhrleistet 
durch die Invarianz des Parameters | 


220 


oe we 


co) 


| : 
= const’. (23) 


In den Vergleichsfallen miissen demnach die Querschnitte proportional zu 
(1 — M2) sein. Unter dieser Voraussetzung sind 


~ 


u(x,7) und v*r7(x,7) 


in den Vergleichsfallen dieselben Funktionen. Fiir Rotationskérper wachst 
in Punkten der Achse tiber alle Grenzen, 1m Gegensatz.zur Profiltheorie, wo 
u(x; 0) endlich bleibt und dem Wert auf der Kérperkontur gleichgesetzt werden 
kann. Daher ist es sinnlos, die logarithmisch-unendlichen u-Komponenten auf 
der Kérperachse zu vergleichen. Deshalb wird hier auf den endlichen Raum- 
einfluss zuriickgegriffen, der in Kérpernahe konstant gleich seinem Wert auf 
der Achse ist. Damit gilt an Rotationskérpern: 


Qz2,(0) OO) Uoo 2 i Qx~(0) 


Whee Qeal*) Inv = f(x) = 2% U(x, y, 2) Sauiico 1 Qea(#) In (7 Qn 2(9) ) ] 


In den Vergleichsfallen sind die auf der rechten Gleichungsseite stehende 
Ausdriicke Funktion von x allein. Zu den Ahnlichkeitsgesetzen fiir Schallnahe 
an Rotationskérpern kommt man mit Gleichung (24), indem man fiir 7 den 
K6rperradius einsetzt, 

Hier sei darauf hingewiesen, dass im Buch Gasdynamik von K. OSwaTITSCH 
(Springer Verlag, Berlin 1952), in den schallnahen Ahnlichkeitsgesetzen auf 
Seite 355 vier gleichartige Druckfehler stehen. Die zweiten Ableitungen des 
Querschnittes — dort F"(x) und K’(x) — miissen umgekehrte Vorzeichen 
erhalten. 

Die Ahnlichkeitsgesetze fiir Schallndhe und Schallanblasung an Fliigeln 
kleiner Spannweite folgen nunmehr aus der Anwendung des Aquivalenzsatzes 
(11) auf das Gesetz (24) fiir Rotationskérper. Nach Division durch Q,,,,(0)/2a 
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erhalt man die analytische Form: 
F) +9(x) Pgh a : 
Fy _f 9% #) InV'y? + (t= 2)? dt] 


—o(*) 


id ) z (C ; 
<a mtn In si = Funktion von x, (25) 


22. [A Y, 2) — Ugo 1 
Qx (9) Uoo 22 


_ oder in Worten: Bei der Umrechnung auf einen Korper affin veranderten Quer- 
schnittverlaufes andert sich die Mach-Zahl der Anstrémung entsprechend 
Gleichung (23) und die w-Komponente entsprechend Gleichung (25) so, dass 
die linke Gleichungsseite die gleiche Funktion von x bleibt. Dabei kann die 
Querschnittform des Fliigels beliebig gedindert werden, wenn nur der Affin- 
verzerrung der Flache Q(x) geniigt wird. 

Von besonderem Interesse ist der Spezialfall einer Affinverzerrung von 
Grundriss und Langsschnitt oder einer der beiden Schnitte. In diesem Falle 
ist Q,,,(0)/2 2 proportional dem Produkt aus Halbspannweite 6 und Dicken- 
verhaltnis 6, was allein schon aus den geometrischen Bedeutungen der Gréssen 
folgt. Aus Gleichung (23) geht dann hervor, dass fiir solche affinverzerrte 
Fliigel folgender Ahnlichkeitsparameter derselbe sein muss: 


pte? z 
Te eg (26) 


Bei allen Problemen kleiner Spannweite kommt es im wesentlichen nur 
darauf an, die w-Komponente an einer einzigen Stelle des Querschnittes zu 
kennen, da dann die w-Verteilung im ganzen Kérperquerschnitt durch ein ein- 
faches Integral (die Querschnittstro6mung) gegeben ist. Bei Fliigeln wahlt man 
am geeignetsten die w-Komponente auf der x-Achse u(x, 0, 0), die praktisch 
der u-Komponente im Mittelschnitt auf dem Fliigel gleichkommt. Auf der 
Achse tritt dann in Gleichung (25) der Ausdruck auf, der mit Gleichung (6) 
umgeschrieben wird: 


o(*) 


ono ox | 1% lel ae 


— a(x) 


+1 
x, t/o) t t 
& i ue lous In o(x) + 2 q(x, 0) a(x) i Me = ae ie ] 


1 
aed we =f abe ts In( say) e Conk 


ZAMP VII/4 
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Der Ausdruck in den eckigen Klammern muss nun bei affinverzerrten 
Fliigeln derselbe sein, denn nach Voraussetzung sind sowohl 


a(x), —_g(%, 0) o() 


a) 


als auch das Integral die gleichen Funktionen von x. Dieser Teil bleibt bei der 
Affinverzerrung ungedndert und kann in Gleichung (25) in die Funktion /(x) 
hineingezogen werden. Ferner ist: 


in-2z2() _ In (b 8) + const . 
2 ot 


Nach einfacher Reduktion lautet das Ahmlichkeitsgesetz fiir affinverzerrte 
Fliigel zusammen mit Gleichung (26): 


1 u(x, 0, 0) — Meo Qn a(*) 3/2 81/2) | _ TF : é 
fo. | ne ee In (Coe Gut ) = Funktion von *. (27) 

Das Ahnlichkeitsgesetz fiir Rotationskérper kann in genau derselben Weise 
geschrieben werden, wenn die Spannweite 25 durch die maximale Dicke t 
ersetzt wird und die u-Komponente an der Ké6rperoberflache y = R(x) ge- 
nommen wird: 


u(x, 0,0) > u(x, R,0); 267, . (28) 


also (entsprechend OswaTITSCH-BERNDT [8], Gleichungen (16) und (17)): 


1 [= R, 0) — Yoo a Qe al*) 


oo as fe Inz?| = Funktion von x. (29) 


Diese Ahnlichkeit der Formeln kann zunichst liberraschen, weil beim 
Fliigel w in der Projektionsflache genommen werden kann, wahrend man beim 
Rotationskérper an der Oberflache bleiben muss. 

Beim Fligel gibt es nach Gleichung (27) eine bestimmte Verzerrung von 
Dickenverhaltnis 6 und Halbspannweite 6, namlich 


b?/? 6¥? — const , 


bei welcher das logarithmische Glied fiir Vergleichsfliigel unverandert bleibt 
und folglich in Gleichung (27) weggelassen werden kann. Zusammen mit Glei- 
chung (26) lasst sich diese spezielle Verzerrung auch schreiben: 


b6'* = const: - - Mr te = f(x); | 


By tS M2, =const; 6-3/4 — M? = Const - | 
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Die untere Zeile ist dabei fiir M., = 1 von selbst erfiillt. Die Verzerrung (30) 
ist bekannt. Sie ist jene Form der Ahnlichkeitsgesetze fiir Schallanstr6mung 
oder Schallnahe, wie sie fiir Fliigel bel‘ebiger Spannweite gilt?°), 

Wie schon frither ausgefiihrt!), werden Stésse in der Theorie kleiner Spann- 
weite im allgemeinen nicht richtig umgerechnet. Dies gilt notwendigerweise 
auch fiir die Ahnlichkeitsgesetze und ist auch hier im allgemeinen von unter- 
geordneter Bedeutung. Die Umrechnung des Stosses im Mittelschnitt ist jedoch 
im Falle von Gleichung (30) richtig, sowie die Umrechnung der Stésse in der 
allgemeinen Ahnlichkeitstheorie fiir Schallnihe stets richtig erfolgt. 


3. Widerstand, Aquivalenzsatz und Ahnlichkeitsgesetze 


Die Betrachtungen dieses Abschnittes seien auf Rotationskérper und diinne 
Fliigel kleiner Spannweite beschrankt. Wahrend die Geschwindigkeit auf dem 
Rotationskérper nicht der u-Komponente gleichgesetzt werden darf und die 
Querschnittstrémung stets an der Kérperoberflache genommen werden muss, 
ist es erlaubt, die Geschwindigkeit auf der Oberflache des diinnen Fliigels 
durch die «-Komponente auf der Grundrissflache anzundhern. 

Mit dem Symbol 2,(x) fiir die w-Komponente des Raumeinflusses gelangt 
man fiir die Geschwindigkeitsstérung eines Rotationskérpers zu folgendem 
Ausdruck auf y = R(x) mit Gleichung (10): 


eR) MR) _ 5 (AY = 0,(0) + Peel ine + 5 (SY 


Moo os los Bk Fag ie 
5 ik Darah) 4 QUE) i Me ee OA) (31) 
Oe ear lr al rp 


Fiir den Fliigel hingegen gilt mit Gleichung (7) und Gleichung (5) in der 
hier angewendeten Naherung: 


+ 0(x) 
ACTS Bat 2 Och, ee / h,(x, #) In| — 2| dt. (32) 
tec, Wes 7 B% J 


Nach dem Aquivalenzsatz sind 2,(x) in Gleichung (31) und Gleichung (32) 
gleich, wenn nur dasselbe Q(x) genommen wird. ae 
Fiir das lineare Problem lautet der Raumeinfluss nach F. KEUNE [2}!*) mit 
_ x9 als Koordinate des Kérperendes (Figur 1) 


WM oo <i: 


Q(x) Qea(*) Tes x (% — *) ce di ee ié Qnalé) dé: (33) 


|~— él 


10) Vergleiche [4], Seite 348 ff., Gleichung (33), und Seite 351 ff. - \ 
11) Vergleiche [7], die letzten 3 Abschnitte des 8. Kapitels. 
12) Gleichungen (4.12) und (4.13). 
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Oe ee Te eet oe i Que ee ie 9 Ae 


Beide Formeln sind fiir alle Werte innerhalb des Kérpers (0 < « < x») giiltig, 
wenn die Ableitung der Querschnittflache Q,(0) in der Kérperspitze ver- 
schwindet. 

Dariiber hinaus ist in Unterschall auch noch Q,(%9) =0 gesetzt im Heck des 
Kérpers (Figur 2c). Es wiirde fiir Q,,(%9) + 0 mit [Q(%9) = 0] rund (Figur 2d) 
oder wie ein Halbkérper [Q(x,) > 0] (Figuren 2a, 2b) endigen. Beides fihrt 
zu derart starken Anderungen in wu, dass die Vernachlassigungen und folglich 
auch die Resultate in der Umgebung von x = x, falsch werden. 

Bei Uberschall hingegen, Gleichung (34), macht sich die Form des Kérper- 
endes stromaufwarts nicht bemerkbar. Der K6rper oder Fliigel kann an irgend- 
einer Stelle abgebrochen und in seinem vorderen Teil als neuer K6rper ange- 
sehen werden, wenn er nur tiberall geniigend schlank ist. 

Endet ein Rotationsk6érper rund (Figur 2d), Q,(%)) + 0 und Q(x») = 0, so 
gilt Gleichung (34) soweit der Korper genitigend schlank ist. Endet er bei aus- 
reichender Schlankheit mit einem Knick, Q,(%9) + 0 und Q(x») + O (Figur 20), ~ 
so gilt Gleichung (34) fiir alle x vor dem K6rperende %, (x < %p). 

Unter Widerstand D ist im folgenden stets die Resultierende der Drucke auf 
die Oberflache des nicht angestellten Kérpers verstanden. Bei abgebrochenen 
K6orpern ist nur der «Mantelwiderstand» gemeint und der Widerstandsanteil 
des Bodens Q(x) nicht beriicksichtigt. 

Der Widerstand D ist gegeben!*) fiir den Rotationskérper 


D W 
Pe me ice is 
0 
und fiir den Fliigel 
_ Pacem 
ga eA pe aa 86) 
; 


In (36) ist tiber die Fliigelflache F zu integrieren. Beide Formeln gelten in 


Schallnahe bis einschliesslich zu Gliedern zweiter Ordnung, sind dort also be- 
sonders genau, 


18) Siehe [4], VI. Kapitel, Gleichung (141). 
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Fir Rotationskérper folgt dann mit Gleichung (31) und Gleichung (10) 


Ovo ey Pine —2 / 92,(x) Qe(x) ad 


= =2 | Q(x) Qalx) dx — 3° = Q2(x) In R(x). 


Fur den Widerstand eines Fliigels tritt nach Gleichung (36) mit Gleichung 
(32) ein Integral auf, welches die Querschnittstrémung enthalt und das nach 
einiger im Anhang (4. Abschnitt, I) durchgefiihrten Integration tiber x mit 
Gleichung (5) die einfache Gestalt annimmt: 


— pot 
Tr=| | hale, 2) 5. | [hel ay tal) ara awa, | 


— G(x) | 


ig (x9) 


oe = fal ae ) fa (Xo, #) In| #2 — 2?| dé dz (38) 
2=0 t~0 
2 (9) 
=; /9 a2 | a %q, t) In| — 2| dé dz. | 
z=0 —0(%») 


In diesem Integral tritt nur die ebene Querschnittstrémung in der Hinterkante 
der Breite 2 o(x,), (Figur 6a) auf. Mit einer Umformung von Gleichung (38) mit 
den Gleichungen (5) und (6) folgt unmittelbar aus den Gleichungen (36) und 
(32) die Formel fiir den Widerstand eines Fliigels: 


D - ee: 
puke = ~ 2 | Ql) Oxla) dx — Ze PR) In ol) 
7) : G(X) +0(X) i <j 
a | ae ) [avo t) In ale.) dtdz, (39) 
2=0 t= —o(%9) 


die in ihrem Aufbau derjenigen fiir Rotationskérper (37) sehr ahnelt. 
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In unterkritischer Unterschallstromung und insbesondere im Fall der lineari- 
sierten Unterschallstromung miissen stumpf endende Fliigel und Rotations- 
korper (Figuren 2a und 26) ausgeschlossen bleiben, da fiir sie die Stromung 
in der Umgebung der Hinterkante falsch wiedergegeben wird. Jedoch muss 
das D’Alembertsche Prinzip unter den Voraussetzungen fiir die Querschnitt- 
flachen Q,(0) = Q,(%9) = 0 der Kérper erfiillt sein. Fiir Fliigel sind diese Bedin- 
gungen nur erfiillbar nach Gleichung (6), wenn g(x, 2) entlang der Hinter- 
kante verschwindet. Damit ist aber fiir Fliigel das Integral /; = 0 und in 
Gleichung (39) bleibt das Integral tiber den Raumeinfluss 2,,(x) allein zu unter- 
suchen. Dieses Integral verschwindet nach den Untersuchungen im Anhang 
(4. Abschnitt, II) fiir beliebige Querschnittflachen-Verteilung Q(x) zwischen 
O<*%< x und Q,(0) = Q,(%)) = 0, so dass der Unterschallwiderstand unter 
den gegebenen geometrischen Voraussetzungen (Figur 2c) fiir alle Aquivalenten 
Korper Null wird, wie erforderlich. 

Dabei sei folgendes hervorgehoben. In 4lteren Arbeiten!4) wurde das qua- 
dratische Glied 1/2 (v/u,,)2 der Geschwindigkeit in der Widerstandsberechnung 
fiir besonders schlanke Rotationsk6érper vernachlassigt. Dieses Glied ist jedoch 
fiir Rotationskérper und insbesondere fiir Halbkérper von gleicher Gréssen- 
ordnung wie die Querschnittstrémung. Es gibt zusammen mit diesem Stré- 
mungsanteil in Gleichung (37) ein vollstandiges Differential, und nur dadurch 


bleibt das d’Alembertsche Prinzip im Unterschall erfiillt» Nur in Ausnahme- | 


fallen, wie bei zur Kérpermitte symmetrischen Kérpern, verschwinden sowohl 
der Anteil von 1/2 (v/u.)? als auch derjenige der Querschnittstromung fiir sich 
allein. 

Fiir den Uberschallwiderstand eines Fliigels im Linearisierungsbereich folgt 
mit den Gleichungen (39) und (34): 


P (Xo) +0(%) 
D _ Qal%o) 2 %q Lf eae) 
me UE /2 a on VM, ni = fa ) a t) A - dt dz 
MeO, ae 
5 al yp aK = | al) g(x) dx , (40) 
0 0 
mit der Funktion 
“oO A 
g(x) = + f Basle — Saal) ge, (41) 


Die in [2], Gleichung (6.15), angegebene Formel muss richtig lauten: 


D 


et [oie ide [O89 SS 


(% — *) 


14) Beispielsweise [11], Gleichung (36). 
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mit einer Funktion G(x), die bis auf die obere Itegrationsgrenze, dort x», mit 
Gleichung (41) iibereinstimmt. Gleichung (40)-geht fiir Q,(x)) = 0 in diese 
Formel iiber und ist im wesentlichen mit der von WARD [1], Gleichung (37) 
angegebenen Formel identisch, wie spater ausgefiihrt wird. 

Ebenso folgt aus den Gleichungen (37) und (34), wieder mit Gleichung (41), 


der Uberschallwiderstand eines Rotationskérpers im Linearisierungsbereich: 


d 


D O24) 2% 


ee eI 
Poo Ug,/2 25 


R(x) (M3, —1 


Da der Bug des Kérpers nach Voraussetzung (Figur 2) wie eine Kegelspitze 
Q(x) ~ x ausgebildet sein soll, ist der Uberschallwiderstand von Fliigel und 
Rotationskérper stets endlich, solange im Heck des Kérpers bei von Null ver- 
schiedenem Q,(%9) auch R(%») und o(%») nicht Null werden), 

Da nach dem Aquivalenzsatz der Raumeinfluss iquivalenter Korper bei 
gleicher Mach-Zahl derselbe ist, folgt ohne Einschrdénkung auf die Linearisierung 
aus den Gleichungen (37) und (39) fiir Uberschall: 


Driiiget Da. R.K. d * R(x) 
Boo Ua /2 Ge u3,/2 ‘ 2m Qz(%o) In (%9) 
1 G(X) +6(%») | F l 
pi ges J Gt de, aS 
= | 12) [ alt) Ino (43) 
s=0 t= —a(%) 
oder 
. a(%») o(%9) R(%) 
Qoo 2 i 
Drags = Yaga. + —Z7 Moo ics 2) | ax t) In jt— z| didz. (44) 
eo t= —o(%,) 


Im allgemeinen ist der Fliigelwiderstand nicht gleich demjenigen seines 
aquivalenten Rotationskérpers. Fiir den Sonderfall ¢(%9, z) = 0 sind beide 
Widerstande gleich. 

Die Widerstandsformeln (43) und (44) kénnen auch aus der bei Warp [1], 
Gleichung (37), stehenden Formel abgeleitet werden, wenn das Potential an der 
-Hinterkante und seine Ableitung in Normalenrichtung bei WARD durch die 
—Quellbelegung ausgedriickt wird. Fiir den Fliigel geschieht dies etwa durch 
unsere Gleichungen (3) und (7). Jedoch ergibt sich bei WaRD das Integral an 
‘der Hinterkante x = x) unter Beschrankung auf lineare Uberschallgleichungen 
aus dem Impulssatz nach einer Integration tiber die Ebene x = x9. Hier ist die 


15) Ausserdem bestatigt diese Gleichung die fritheren Ausfilhrungen, wonach man zu einem 
fehlerhaften (weil unendlichen) Widerstand kommt, wenn R(x9) = 0 wird bei Qy(%9) + O (Figur 2d). 
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Integration mit Hilfe der Geschwindigkeitsverteilung tiber den ganzen K6rper 
ausgefiihrt worden. Die Gleichungen (43) und (44) gelten deshalb bei schall- 
nahen Strémungen auch fiir abgebrochene Kérper (Figur 2b), wenn das 
Korperende geniigend tief in Uberschallstrémung steckt, so dass von dort keine 
Stérungen mehr stromaufwarts wandern kénnen. Ausgenommen bleibt stets 
der Kérper vom Typus Figur 2d. 

Der Widerstand aquivalenter Kérper unterscheidet sich nach Gleichung (44) 
nur durch die Querschnittformen und die aus ihnen folgenden Querschnitt- 
strémungen im Endquerschnitt des Kérpers. Dieser Unterschied ist unabhangig 
von der Mach-Zahl, so dass daraus eine niitzliche Formel folgt, in welcher der 
Index L eine beliebige Mach-Zahl-Linearisierung in Uberschall kennzeichnet: 


hice Paks Drriiget FF Dig RK. = Dy rriigel ae Dy ag RK. . (45) 


Die Differenz der Widerstande aquivalenter Koérper ist von der Mach-Zahl un- 
abhangig. 
Gleichung (45) ist trivial fiir unterkritische Strémung, da dort jeder Sum- 
mand der Gleichung schon wegen des d’Alembertschen Prinzipes fiir sich ver- 
schwindet. Im Uberschall kann die Differenz der Widerstande jedoch positiv 
oder negativ sein; dann handelt es sich um KG6rper (etwa Figur 20), fiir welche 


die Unterschallvoraussetzungen verletzt sind. Gleichung (45) verliert ihre © 


Giltigkeit, sobald ein solches Kérperende in Unterschallnahe gerat. 

Die Gleichungen (43) oder (44) und (45) geben nicht nur den Gesamtwider- 
stand aquivalenter Kérper, sondern bei ausreichender Schlankheit den Tezi- 
wtderstand von der Kérpernase bis zu einer beliebigen Stelle x (x > x) inner- 
halb des Kérpers. Dann sind aber im allgemeinen die Anderung Q,(x) der 
Querschnittflache und die Querschnittstrémung von Null verschieden. Die 
Teilwiderstande aquivalenter Kérper unterscheiden sich — wie die Geschwindig- 
keiten, Gleichung (11) — durch die Querschnittverteilungen im betrachteten 
Querschnitt. Selbst bei gleichem Gesamtwiderstand von Fliigel und aquiva- 
lentem Rotationskérper, also fiir Q(%»9) = Q,(%9) = 0, sind im allgemeinen die 
ortlichen Widerstandsanteile des Fliigels und des aquivalenten Rotationskér- 
pers verschieden. Gleich sind sie nur in Ausnahmefillen fiir den Teilwiderstand 
bis zur Stelle des Maximalquerschnittes, in dem Q,(x) = 0 und Q(x) > 0 sind; 
wenn dort die Quellverteilung (x, z) tiber die értliche Spannweite Null ist, die 
im allgemeinen positive und negative Werte annehmen kann. 

Mit den vorliegenden Ergebnissen sind nun auch die schallnahen Ahnlich- 
Reitsgesetze fiir den Widerstand gewonnen. Mit dem auf die Fliigeltiefe als Lan- 


geneinheit bezogenen Maximalquerschnitt Qmae folgt fiir Rotationskérper nach 
[8], Gleichung (18b): 


PRK Roam, te ai Qz(%o) 
Or att wn 2 Yoo 29 Q? ax MO nee = const . (46) 
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Hier kommt man zu dem gleichen Ergebnis, wenn Gleichung (24) in Glei- 
chung (35) eingefiihrt wird und beachtet wird, dass in W/u.—1 das Glied 
v?/2 nach Gleichung (20) proportional Q,,,(0)/2 a ist. Schliesslich folgt Glei- 
chung (46), da Q,,(0)/22 ~ Q2,,, ist: 

Der Quotient Dp x./Q%,., des Widerstandes im Ahnlichkeitsgesetz deutet 
direkt auf den in der Ballistik iiblichen mit der Stirnfldche Q naz dimensionslos 
gemachten Widerstandsbeiwert hin: 


a) 2 
D= a Wo ae Ca- 


Fir aquivalente Koérper mit affinverzerrtem Querschnittverlauf ist mit 
Gleichung (43) o(x9) proportional der auf die Gesamttiefe bezogenen maxi- 
malen Halbspannweite b und R(x») proportional VQ pax: 


ee TS 2 2 
Drriiget a Dig rx. rT Sco. Un 2 x Q(X) in /2nae + const Onan , (47) 


oder 
Qo ,2 1 Az, 
Drtisget a ne Uno Qn Q2(%) Ind 
Qo 2 1 ap /O- 2 
ay Dig rk. * ie Ux os) Q:(%o) In) ads -- const Oe. : 


Damit nun auf der rechten Gleichungsseite die schallnahe Ahnlichkeit, Glei- 
chung (46), steht, ist zunachst auf beiden Seiten ein Ghed mit InQ,,,, zu 
addieren und durch Q3,,, zu dividieren: 


DFiiiget Oc 2 1 (92 (%.) ae 
8 a = 2 oo 22 a. O2 ax Ind VO cic 

Dig.R.K. oa ae Qz(%o) at Pi a once 

—y, Ofnax 5 2 Yoo 22 Oz ax Ona =} ( ) 


Vergleicht man bei gleicher Mach-Zahl nur Fliigel miteinander, so ist mit 
Gleichung (46) die ganze rechte Gleichungsseite eine Konstante und 


Onan ~ bu 6: ae + 22 12, — oe. In b32 6/2 = const. (49) 

Dieses Ergebnis folgt natiirlich auch, wenn in Gleichung (36) das Ahnlich- 
keitsgesetz fiir Fliigel, Gleichung (27), eingesetzt und die Integration ausge- 
fiihrt wird. Mit b ~VOmax geht Gleichung (48) tiber in Gleichung (46), ebenso 
wie die schallnahen Ahnlichkeitsgesetze fiir Fligel, Gleichung (27), in die fir 
Rotationskérper, Gleichung (28), tibergegangen sind. 
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Beim Ubergang von einer Mach-Zahl auf die andere gilt Gleichung (26) 
fir Fliigel und mit 6b ~ Qa» fiir Rotationskorper, das heisst, 1 — M2, ist 
proportional zu Qjyqx Zu andern. 

Die Widerstandsformeln seien an folgenden Beispielen erldutert. Zuerst sei 
der Fall betrachtet, dass die Quellen im Endquerschnitt verschwinden 
[¢(%», 2) = 0]. Dann ist nach Gleichung (6) auch Q,(%9) = 0, der a4quivalente 
Rotationskérper endigt entweder als Halbkérper wie Figur 2a oder ais Voll- 
korper kegelig wie Figur 2c. Die Aquivalenten Fliigel enden dann ohne Knick 
in der Hinterkante, bei Vollfliigeln also wie die Langsschnitte in Figur 3. Die 
Theorie kleiner Spannweiten ist in diesem Fall nicht nur in Uberschall-, sondern 
auch in Unterschallstrémung tiber die Hinterkante hinaus giiltig. Der Wider- 
stand ist in Unterschall Null und in Uberschall nach Gleichung (43) fiir alle 
Aquivalenten Korper gleich und nach Gleichung (47) bei Affinverzerrung pro- 
portional zu Q%,,,. In den Teilwiderstanden unterscheiden sich solche Korper 
allerdings. 

Diesen sei ein «abgebrochener» Fligel nach Figur 5 gegeniibergestellt, dem 
ein aquivalenter Rotationskérper vom Typus Figur 2b entspricht. Die Unter- 
schalltheorie kleinen Seitenverhaltnisses ist in naherer Umgebung der Hinter- 
kanten hier nicht mehr anwendbar. Die Quellverteilung in der Hinterkante 
zeigt im Mittelsttick Senken, an den ausseren Randern jedoch Quellen ent- 
sprechend der Dickenabnahme und -zunahme an diesen Stellen. Ferner ist 
Q.(%) + 0. Damit unterscheiden sich die Widerstande von Fliigel und aqui- 
valentem Rotationskérper entsprechend Gleichung (43) oder (44); bei Affin- 
verzerrung entsprechend Gleichung (47). Auch die Teilwiderstande unterschei- 
den sich im allgemeinen selbst im Punkt des Querschnittmaximums Q,(x) = 0. 
Die Geschwindigkeitsverteilung im Linearisierungsgebiet von M., > 1 ist in 
[3], Bild 10 bis 12, gegeben. 


4. Mathematischer Anhang 


I. Der Beweis fiir Gleichung (38) soll in méglichst allgemeiner Form durch- 
gefiihrt werden, indem 


a) an Stelle der Dickenadnderung h,(x, z) nach Gleichung (5) die Quellver- 
leitung 9(x, z)/2 geschrieben und 
b) fiir diese und den Grundriss des Fliigels (Figur 6a) nur Symmetrie zur 


x-Achse, q(x, +2) = q(x, —2), vorausgesetzt wird. Dann kann die Integration 
uber ¢ auf eine Fliigelhalfte beschrankt bleiben: 


a(x) 
th=/ [axa 2 Fate In| t? — 22| a dj. 
F 


0 
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Die Funktion In | ¢? — z?| ist eine in ¢ und x symimetrische Kernfunktion 
fo(t) 2) = K(z, t), ferner konnen nach Figur 6a fiir die halbe Belegungsbreite 
2. = s(x) die Umkehrfunktionen x = ¥,(z) und x = x4(z) eingefiihrt werden, fiir 
welche gilt: 


O0<2< db: z=o(*,(2)); und a<z<d: z= o{ x,(2)) . (50) 


y/ 
/s _A[x,z =konst) 


Alxz)=0 
Figur 5 Figur 6 
Deltafliigel mit endlicher Dickenverteilung Fligelgrundriss (a) und Skizzen 
in der Hinterkante. Die Dickenverteilung ist zur Umkehrung der Integrationsgrenzen. 


durch geradlinige Erzeugende konstruiert und 
an den seitlichen Fligelrandern Null. 


Dann ist /; in grésstmoéglicher Allgemeinheit, iiber eine Fliigelhalfte integriert : 


Oy Xo a(x) 
Bh [ae 35 fa. t) K(t, 2) dt} dx dz 
2=0 *=%,(2) 0 
b tals) a(x) 
+ ¢ g(x, 2) 2 p g(x, t) K(t, 2) dt| dx dz, (51) 
0 


Oy x= %4(2) 


und es ist zu beweisen, dass mit der Belegungsbreite 6) = o(%») 


6(%o) o(%s) 
2n-4 [ [Kea ase 2) qr) dide=0=2I,— 2h, (2 
z=0 ¢=0 
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ist, wenn die erste Form von Gleichung (38) gelten soll. Der Beweis wird hier so 
skizziert, dass seine Nachpriifung ohne Umwege durchgefiihrt werden kann. 
Wird zunichst Gleichung (51) partiell nach x integriert, so ergeben sich drei 
Doppelintegrale und zwei dreifache Integrale. Von dem einen Doppelintegral, 
namlich 


dy 
as 


Df ee / | K(t, 2) q(%q, 2) q(%, t) at dz + usw. , 


werden soll. Deshalb ist die Hilfte dieses Wertes auf die linke Gleichungsseite 
von Gleichung (52) gezogen, und in den drei Doppelintegralen J}; nach der 
partiellen Integration bleibt die andere Halfte stehen: 


% % b z 
2 I= [Kt2) avez) geet) atde + | | Kb) 9(xale),2) a(xal2),#) at da 
Ole) Z=6,t=0 : 

= 
—/ [Ke 2) q(44(2), 2) @( (2), #) dt dz. 


Z—017=0 
Fiir die beiden dreifachen Integrale nach der partiellen Integration von (51) gilt: 


Cpe 1£ a a(x) 


a 


Cig [ ic 2) | q(x, t) K(t, 2) dt dx dz 


Z 50 (2) 0 
b  x2(2) a(x) 


+/ ri qu(%, 2) / g(x, t) K(t, 2) dt dx dz. 


2=09 *,(2) 0 


Diese beiden dreifachen Integrale werden umgeformt, indem zundchst die Intal 
grationsreihenfolgen umgekehrt werden: 
1. nach Figur 66 die von ¢ und x: 


a) MO soe Ss 


Xp a(x) Xo Xs b x9 (t) 


Jets) atax =| Jet mavars [ [ss 9 avats | [sis axat, 


x,(z) t=0 t=0 x,(z) t=2 2,(t) t= %,(t) 


(2) o(x) Z  #%9(2) b  x9(t) 


. 
2=07t=0 
stimmt der halbe Wert mit Gleichung (38) tiberein, wahrend der Rest zu Null 


bio, S28: | [etx dtdx— [el 9 avaes [ [axdt. 
%,(z) t=0 


t=0 x,(z) t=z x,(t) : 
2 h Fi i ti 
- nach Figur 6c in den nun fiir 2 J; entstehenden fiinf Integralen die von : 
t und z, mit folgender Vertauschung der beiden Integrationsvariablen z und 4: 
: 
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a) im zweiten Integral: 


Sy % ai “ae 
Jd g(t, 2) dt dz alg g(t, 2) de dt of fat t) dt dz, 


b) im dritten Integral: 


o b Gy 
it “feu 2) dt dz = -f / ste 2) dz dt > [ [se t) dt dz, 
2=0 a, =o, 2-0 z=) t=0 

c) im letzten Integral: 

b ob ee: b 2 
/ eG 2) dt dz =| 1G z) dz dt >/ Jee t) dt dz, 
Z=0) 2 t=, 2=06, Z=0) t=0 


wahrend das erste und vierte Integral unverandert bleiben. Nun sind im 
dritten und letzten der fiinf Integrale die beiden Integrationsbereiche tiber ¢ 
bei gleichen Integrationsbereichen z und x zusammenzufassen zum Integra- 
tionsbereich 0 < ¢ < z des dritten Integrals mit dem Ergebnis in unveranderter 
Reihenfolge: 


Go x; 


2Iy.= | [Keo t) | (au( x, 2) q(x, t) + qo(x, t) q(x, 2)] dx dt dz 
z=0 t= “o x,(2) 
5 Zz %2(2) 
+/ | Ke, t) | [du(x, #) g(x, 2) + delx, 2) a(x, t)] dx dt dz. 
<0, to %,(z) 


Da nun aber in eckigen Klammern d/dx [q(x, 4) g(x, z)] steht und die Kern- 
funktion K(x, t) = K(t, z) ist, folgt, nach der Integration tiber x, Zusammen- 
fassung der beiden Integrationsbereiche tiber z der Integrale mit der unteren 
Grenze x,(z) und schliesslich Vergleich mit 2 J;, nach Gleichung (52): 


2 Jn —2 he 


z 


-3/ | K(t, 2) q(% 2) q(%» #) dt — / K(t, 2) ¢(%q 2) q(%q t) di! de. 


t=0 
Nun sind aber diese beiden Doppelintegrale gleich), so dass unmittelbar folgt: 


2hi1—2h2=9, 


16) Vergleiche fiir das erste Integral die Umkehr der Integrationsreihenfolge und Vertauschung 
der Variablen in 2. (a). 


,- 
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womit Gleichung (52) bewiesen ist und im Fall von Gleichung (38) gilt: 


G(Xy) O(Xy) 
‘ 2 
R= 4 [ fale) alee t) inlet — 2] at de 
z=0 t=0 
(4%) +(%9) 
Ze eet (Sere ) [at ) In|t — 2| dt dz (53) 
t = —a(X) 


oder mit der Umformung 2. (a) in mathematisch am besten ausfiihrbarer Form: 
o(%9) Zz 
= y [a2 [ a% #) Ine? — #) dt de. (54) 
2=0 t=0 | 
II. Fir den Beweis, dass der Widerstand, Gleichung (39), im Fall der 
linearisierten Unterschallstromung verschwindet, wird der Raumeinfluss des 
Potentials am besten in der Form geschrieben?’) : 


Q(x) = ri | Quit ) n(x — 8) d+ 1 | Ould Ain 2 (£ — x) des 


wobei fiir diese Formel und die folgende Ableitung Q,(0) = Q,(x)) = 0 voraus- 
zusetzen ist (Figur 2c). Dann lautet das aie Gleichung (39): 


Tin = | Qal2) Qul) dx = — je ) Qoal 
=a, / Li ) Qualé) In2 (x — &) dé ax 


ier f [Qe ) Orxte) &) In2 (€ mad x) dé ax = 


Im ersten Doweinteln ae nun die Integrationsvariablen x und € zu ver- 
tauschen, wahrend im zweiten Doppelintegral die Integrationsreihenfolge ver- 


tauscht werden soll, was mit x = z, x, = 6, é=¢ ganz analog zu Figur 6c zu 
geschehen hat. Damit folgt: 


fume f [Oat ) Qea(x) In2 (E ~ x) dx dé 


§=0 x= 


=e f fdas ) Oral) In2 (& — x) dx dé =0, 


€=0 x= 


fiir beliebige Q(x)-verteilung, wenn Q,(0) = Q,(x,) = 0. 


™) Vergleiche [2], Gleichung (4.11). 
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Summary 


In transonic flow the velocity disturbance on low aspect ratio wings can be 
reduced to the velocity disturbance on bodies of revolution with the same 
distribution of cross-sectional area. This law of equivalence has been proved in 
another paper and is stated here from the physical point of view. By the similarity 
law for bodies of revolution corresponding similarity taws for low aspect ratio 
wings are given. Further the law of equivalence for the pressure drag is given 
for bodies of revolution and pointed low aspect ratio wings with any cross-section 
at the end, if the end lies sufficiently far in the supersonic flow. The last one is a 
generalisation for transonic flow ofa formula of WARD originally valid for linearized 
supersonic equations only. 


(Eingegangen: 16. Dezember 1954.) 
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Zur Auflosung linearer homogener Differentialgleichungen 
2. Ordnung 


Von JosEF ZBORNIK, Chur?) 


Die Differentialgleichung 
gt G pase tee (1) 


mit den Koeffizienten G = G(x) und H = H(x), welche in dem betrachteten 
Intervall stetig und # 0 sein sollen, geht durch die Transformation 


u= fexp(—[Gdx) dx*), ou) = y(x) (2) 
in die reduzierte Differentialgleichung 
u** (v')? + Hous 0 (3 


iiber, wobei *, **,... die Ableitungen nach w anzeigen und H und w’ mittels 
(2) als Funktionen von w dargestellt sein sollen. 
Ist v = F(u) das allgemeine Integral der Differentialgleichung (3), so stellt 
gemass (2) 
y = F({ exp(—/ G dx) dx) (4) 


die allgemeine Losung der gegebenen Differentialgleichung (1) dar. 

Im weiteren werden aus der beliebigen Zahl von Differentialgleichungen, 
die nach dem angefiihrten Verfahren lésbar sind, einige Typen herausgegriffen 
und jeweils die gegebene Differentialgleichung (1), w und die allgemeine Lésun 
(4) angefiihrt?). 


1) Calanda S. A., Domat-Ems. 

*) Die Integrationskonstanten sind bedeutungslos, weil sie in der allgemeinen Lésung de 
jeweiligen Differentialgleichung in den Konstanten C, und C, enthalten sind. 

3) Woeinzelne Differentialgleichungen unter den Beispielen des Werkes von Prof. Dr. E. KAMKE 
Differentialgleichungen, Losungsmethoden und Lésungen, Bd. 1, 4. Aufl. (Akademische Verlags: 


gesellschaft, Leipzig 1951) auftreten, wird durch Angabe der betreffenden Beispielnummer darau 
hingewiesen. 
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1. Differentialgleichungen mit der reduzierten Differentialgleichung 
v**+qu=0 (q= const) 


v= Ci exp(i wg") + Cy exp (—é w gt”) = c, sin (wg!) + cy cos(w q¥?) 


y" = (5) y' — B27? y=0, f= f(x), {+0 [Kamxe (2.79)]. (1.1) 


w=f, y=Cyexp(bf) + Cyexp(—0f) = cy sin (ib f) + cy cos (i bf). 


y"—-(F + a£) y'— vg py—o [f= Hx) +0, g=g(x) +0]. (1.2) 


g 


m= [ ferdx, y = C, exp (bu) + C, exp(—b u) = ¢, sin(i b u) +c, cos(i bu). 


y" — cy’ — b% ey —0 [Kamxe (2.33, 2.34)], (2.3) 


u—e*/c, y=C,exp(bu) + C,exp(—bu) =¢, sin(t bu) +c,cos (tb u). 


" / J f , es . , 
y + (ate E\y aq y= [f= f(x) = +i, f +0}. (1.4) 


w—In(f+ (P+ 4%)""), y= (f+ (P+ a)?)? + (f+ (+077) 


= ¢,sin(? bu) + cgcos(t bu) . 


Fir 
Ge? b= 5 (1.44) 


erhalt man die Differentialgleichung [KAMKE (2.81)| mit der Lésung 


y= a (f de (f? s% RP ue ae Cc. (f Ee (7? a 4)1/) -1/3— 


Wird speziell C, = 1: |2"3| und C, = —|2”?| angenommen, so gilt y* + 3 y = f. 


(x?+2Ax+ B)y"+(x+A)y’-—my=0. elie) 


u=In(x+A+ (x2+2Ax-+4B)”), y= C, (xt At (x? +24 %+ BI? )m 
D4Ca (e+ A+ (x?+2A%+4B)")-™. 

Fir m=1: (2n+1), »=0,1,... und CQ, =C, (A? — B)™ gentigt y der 
- Gleichung 

yem+tt ay Cy Cy y?"— 2 + ay (Cy Caater Pare ae aatGu Caley 


=2C}"+1(x+ A), 


ZAMP VII/5 


66 JosEer ZBORNIK ZAMP 


wobei die Koeffizienten a, (k = 1, 2, ... , m) durch die Rekursionsformeln 
om ‘2n+1 2n—1 2n—3 ‘2n—2k+3 
vem = 2) 9a) <n 


bestimmt sind. 


(27 x2 +4) y" + 27% y'—-3y=0 [KamKE (2.290)], (1.6) 
Sonderfall von (1.5) mit A=0, B=4/27, m=1/3; a,=-—3: 


Nyse C, E ae (x9 + ees aed pare E ra (22 4 ee 


Fir Cy=—1:|2"*|, C,= 1/3 |2**| gilt 42+ y+2=0. 


50 x (%—1) y"4+25(2%—1) y’—2y=0 [Kamge (2.292)|, (1.7) 
Sonderfall von (1.5) mit A=—1/2, B=0, m=1/5, —a@=4@=5: 


1/5 —1/5 


y=Cy (« = es + (x? — xy") Ere (x es = + (x%— xy) 


Fir C,=1/C, = |2**| gilt y8—5y9 +5 ye4eo2. 


(x?+ 1) y"+a%y'+2y=0 [Kamne (2.222)], (1.8) 3 
Sonderfall von (1.5) mit A=0, B=1, 
y=C,sinw+C,cosw, w= 2? In(x + (x?4.1)'). 


(x? — x)? y” + (2% — 1) (x? — x) y’ —my=0, (1.9) 
u=In((x—1)/x), y= C,exp(mu) + C,exp(—mu) = C, ((x — 1)/x)™ 
+ Cy (%/(% —1))™. 


2 SU xX CH XK 


gb RP es yt n* yy dn? =O | KAMRE (2.74) eee 

u= [dn x dx, y=csin(n | dn x dx) + c.cos(n [ dn x dx) . 
gS TCP OAG CE i aR ah tae (1.11) 
u=fexpf dx, y = C,exp(t bu) + C,exp(—i bu) = c, sin(b “) + cy cos(b u). 
y"—metex y+ b2sin?™x y= 0, (112) 


Sonderfall von (1.11) mit f= mInsinx, u—= f sin™ 4 dai, vy 


= ¢, sin (bu 
+ ¢, cos (bu) . moe 
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y" —ctgx y'+sin?xy=0 [KaMKE (2.69)], (1.13) 


Sonderfall von (1.12) mit m=b=1, y=, sincosx + c, coscosx. 


y+ mtgxy’ + b2 cos?™x y= 0, (1.14) 
Sonderfall von (1.11) mit f= mlncosx, u= f cos mde, | =e, Sin (6 4) 


+ c,cos(b u) . 


y" + tex y' — cos*x y=0 [KAmKE (2.67)]. (45) 
sonderfall von (1.14) mitm=1, b=7, y=C,e&"*%+C,e7%". 


sina xy" — 2m y+ 2 b%tg?™-1esinin y= 0, (1. 16) 
Sonderfall von (1.11) mit = mIntgx, u= [tg™wxdx, 
vot, sin (b [ tg™ x dx) + c,cos(b : te" dx): 


sin2 xy" — y'+ 2 b%sin?x y=0, (1.17) 
Sonderfall von (1.16) mit m=1/2. w=/tg’?xdx, y=c,sin(b / te’? x dx) 
+ ¢,cos(b [ te”? x dx) . 


xing y= 4 — a? x ln?" <x y= 0, (1.18) 
Sonderfall von (1.11) mit f=In(In"x), b=7a, 
u= [In"x as = a (—1)"-* wi/klIn®¥x, y= C,exp(au)+ C,exp(—a yu). 
k=0 


xInx y"” — y'’—xIn?xy=0 [Kamxe (2.412)], (1.19) 


Sonderfall von (1.18) mita=n=1, y=Cy, (x/e)”+ Cy (e/x)*. 


xt y" 4+ 2x3 y'4+a%y=0 [KamKE (2.350)], (1.20) 
w= —1/x, y=c,sin(a/x) + ccos(a/x) . 
y" —nIgx y’ — a? Wof?*xy=0, (1.21) 


Sonderfall von (1.11) mit f= InGofx, w= [Cojfrxdx, 
y= C,exp(a [ Cof"x dx) + C, exp (—a f Gof" x dx) . 
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y" = 2cteZ # y! —sin?2 + y=0, (1.22) 


u=1/2sin?x, y= C,exp(sin?z) + C, exp (— sin 3aa)i 


(x2 —1) y"4+4y'-ay=0 [Kamxe (2.235)], (1,29 


Spezialfall von (1.5) mit d4=0, B=—1, m=a@ 
y = Cy(% + (x? — I)" + C(x + (2 — YP). 


02" hy" +E (n+ 1) art nal ay —vW eel" y=0, (1.24) 
u = —1/an exp(a/x"), y = Cy exp(e*’*”) + C, exp(—e7™) . 


y" — (a+ 2cx)y'—Bexp(2ax+2cx*)y=0, (1.25) 


u=exp(ax+cx*), 
y = C,exp(b exp(a x + ¢ x*)) + C,exp(—b exp(a x +c x). 
y" — (F424) y'- (ve f2- a4 —a*) y=0 [Kame (2.80)], (1.26) 
f=f(x), f +0. Die Differentialgleichung geht durch die Transformation — 
yie=e** 2 in (1:1) mit.z-statt.y tiber. y = e** [C, exp(O fi. Cyemp (pane 
ini ae Ge , we ee! — fA o2@ 72 2. salt 
y ie +a£ +20) y+ |e (4 +a) b* g P+ely=0, (27) 
f= f(x) +90, g=g(x) +0. Die Differentialgleichung geht durch die Trans- 
formation y = e°* zin (1. 2) mit z statt y tiber. 
y= ¢*[C,exp(b fg" fax) + Cyexp(—b f g* fax). 


62 yl" — (x5 +20) xy! (—syrar+ct x4 c2 +c) y= 0... (1-28) 
Die Differentialgleichung geht durch die Transformation y = x¢ z in (1.1) mit 
z statt y tiber. 

y = #[C, exp(b f) + Cyexp(— 6 f)] = ¥°[o, sin(é bf) + cy cos (ib f)]. 
Kr yh — [x (4 +24) +26] ay 
ae [22 (a+ at bei) + x(2ac+ c f, 


(1.29) 


)+e(e+H]y=0, 
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f=f(x),f +0. Die Differentialgleichung geht durch die Transformation 

y = x°e" zin (1.1) mit z statt y tiber. 

y = x° e#(C, exp (b f) + C, exp(—b /)]= x et*[e, sin (i bf) + cy cos(i bf)]. 

2. Differentialgleichungen mit der reduzierten Differentialgleichung 
v** w—n(n-1)v=0 


v= Cu" + C,u'-™ firn + 1/2; v=u'? (C,+ C,Inw) firn=1/2. 


y"— (5) 9! n (n (FY y=0, (2.1) 
ens) = O05 f+ 0, w=}. 
Y= C+ Gf, n#1/2; y=P2(C,+Gnf, »=1/2. 


y" +tgx y’ —n(n—1)ctg*xy=0, (2.2) 
Sonderfall von (2.1) mit f= sinx =w. 


y = C,sin®« + C,sint!-"x, »+1/2; y=sin"?x (C,+ Cylnsinx), »=1/2. 


sin?2 x y" — 4sin?xsin2 x y'—4n(n—1)y=0, (253), 
Sonderfall von (2.1) mit f=tgx=u. 
y= C, teg™~+Cytgi-"x, m+1/2; y=tg’?®x(Cy+C,Intgx), n=1/2. 


xy" —(2x2+4+1)y'—4n(n—-1) #8 y=0. (2.4) 
uw = 1/2 exp (x?) . 


B= © e** + C,c8-9*, naij2; ys e2"* (4,7), n= 1/2. 
wT ie (ae = iS 2 in =0 2.5 
fae + 2a) y [7m (n 1) P? a atl y G ( ) 
Die Differentialgleichung geht durch die Transformation y = e**z in (2.1) mit 


zstatt y tiber. 
y= et (Cyt CoP), m+1f2; y=er fi? (C+ CInf), n= 1/2. 


ety” — (x C + 2c) xy'+|—n (w—1) ts tet xte(o+ily=0. (2.6) 
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Die Differentialgleichung geht durch die Transformation y = x°z in (2.1) mit 


z statt y tiber. 
y= xe (Cy fr+C, fir), nelf2; y= xe fl? (Cy + Cylnf), n=1/2. 


wey" — [x (+24) +2¢] cM | 
-} | (a+ f n (n —1) r) +x(2ac+o4,) (2.7) 


+e(c+1)] y=0. 


Die Differentialgleichung geht durch die Transformation y= x? e*”z in (2.1) 


mit z statt y tiber. 
y = x 8 (Cif? + Cy ft-"), n+1/2; y= x erp? (C,+C,lnf), »=1/2. 


x2 y" + xy’ —n(n—1)In-*xy=0. (2.8) 
Sonderfall von (2.1) mit f=Inx=u. 


y=C,In"x+C,Int-"%, n+1/2; y=In"?x(C,+-Cylminx), n=1/2.) 


y" — (L+e) y'—n(n— 1lhe*y=0. (2.9) 
Sonderfall von (2.1) mit / = exp(e) =u. 
y = C,exp(ne*) + Czexp((1— 1) e*), n+1/2; 
y=exp (1/22) (C,+C,2), n=1/2. 


y" — (cosx — tgx) y’ — n(n —1) cos?x y=0. (2.10) 
Sonderfall von (2.1) mit / = exp(sinx) =. 
y = C, exp(n sin x) + Czexp((1— 1) sinx), n+ 1/2; 
y = exp(1/2 sinx) (C,+ C,sinx), n=1/2. 


i] ge , , t : 
y'— (2 +e) yn -1 gr y=0. (2.11) 
Sonderfall von (2.1) mit f=expg=u, g=g(x), fie On 

y = Cy, exp(n g) + Cy exp((1 — n) g), n+1/2; 


y =exp (1/2 g) (C1 + C,g), n=1/2. 
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cos*x y" + cos*x (2 tgx cos?% —1) y= » (n-1) y=0. (2.12) 
Sonderfall von (2.1) mit f = exp (tg x) = w baw. von (2.11) mit g = tex. 
y = Cy exp(ntgx) + Czexp((l— ni) tex), n+ 1/2; 
)y = exp(1/2 tgx) (C,+ Cytgx), n=1/2. 


=. Differentialgleichungen mit der reduzierten Differentialgleichung 
v**+ 6? u?"—2v=0, b= const 


v= ul? Z,., (b/rw), Z= Zylinderfunktion. 


” 


Seay eps pt Ry 0, 6.1) 


f=f(x), PO. w=f> y= ff? Zo, (bir). 
y" — y+ Be y= 0. (3.2) 
Sonderfall von (3.1) mit u = e*, r= —1: y = e*[C,sin(b e~*) + C, cos(b e~*)]. 
y+ y'4+ ett y=0. (3.3) 
Sonderfall von (3.1) mit « = e*, r= —1: y = e*[C, sin (6 e*) + C, cos (6 e*)|. 
y" — ay’ + a2 b? e2tat y= 0. (3.4) 
Sonderfall von (3.1) mit w =e: y=exp(a x/2) Z,)., (b/r e"%”) . 
y" + tex y’ + bd? cos*xsin**?-2x%y=0. (3.5) 
Sonderfall von (3.1) mit w=sinx, y=sin"? x Z,.,, (b/r sin’) . 
y" — Igx y’ + b?Coj?2x Sin?" 2x y=0. (3.6) 
Sonderfall von (3.1) mit w= Ginx, y= Gin"? x Z,., (b/r Sin”) . 
xy" — (142 27) y' +402 x8 ctr y=0. (ee 
-Sonderfall von (3.1) mit «= exp(x?), y = exp(¥?/2) bile ci) 
xy" + xy’ + bln? Ax y=0. (3.8) 


- Sonderfall von (3.1) mit w=Inx, y=In'"? x Z,),, (0/rin"4) . 
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yy" — (f +2 a) y! + (0° f2r2r-24 ata r) yO (3.9) 


Die Differentialgleichung geht durch die Transformation y = e*” z in (3.1) mit 
2 statt y tiber, y= ef” Z,,, O77) - 


wy! — (nF 420) xy’ + (srt prt atte ster +c) y=0. (10) 


Die Differentialgleichung geht durch die Transformation y = x¢ z in (3.1) mit 
zstatt ytiber. y= 2 f" Z, 0) (bir 7"). 


x2 y" — E (F- + 2a) + 2c La 
ies ; a” lee 
2 [72 ue D fe for 2 pls aL 
+ |x (a maar aur ar! )4+x(2actet : (3.14) 
+e(c+1)|=0. 
Die Differentialgleichung geht durch die Transformation y = x° e**z in (3.1) 
iit 2 statt-¥ uber. y= ceo") Zon), 


4. Differentialgleichungen mit der reduzierten Differentialgleichung 
v** (e?2— 1)—4e74v=0 


v = (e2”—1)[C, + C,In(l —e- 241+ C,. 


y ee if ert iA 
yy ty y Peyeat) ae (4.1) 
u=f, y= (ef —1)[C,4+ C,In(l—e-*)]+C,. 
x (%®™ —1) y" + (42™—1) y’ —4n2 x2"-1ly— 0, G2 
Sonderfall von (4.1) mit w= nInx. 
y = (x? —1)[C, + Cyln(l — x-2%)] 4 C,. 
x (x*—1) y" + (x2? -1) yy’ -—4x4y=0. (4.3) 
Sonderfall von (4.2) mitn=1. y= (x?—1) [C,+ C,n(1 — £2) | Ge 
cos*% y" — ctgx y’-—4y=0. (4.4) 
Sonderfall von (4.1) mit w= —Incosx. y = tg" ~ (Cy + 2.C, nsin'#)iaGae 


: 
: 
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(1 — 62%)? y” —2 (L — €2%) y" =F et y 0, (4.5) 

Sonderfall von (4.1) mit « = —1/2 In (1 — e2*) . 

y = exp (2 x)/(1 — e%*) (Cy + 2Cy x) Cy. 


Durch die Transformation e?* > x geht die Differentialgleichung (4.5) iiber in 
%(*—1)? y"+ x (x-1) y'—y=0 (4.6) 
mit der Lésung y = x/(v —1)[C,+C,Inx]—C, [Vgl. Kamxe (2.326)]. 


, 


” : ae 
Sin x Coj x 


y say 0. (4.7) 


Sonderfall von (4.1) mit w=InGofx. y = Gin?x (C,4+ 2C,Tqx) + C,. 


5. Differentialgleichungen mit der reduzierten Differentialgleichung 
v** sin’u — 2 v= 0 bzw. v** cos?u— 2u= 0 


v= ctgu (C,— C,u)+C, bzw. v=tgu(C,+ Cau) +C,. 
ex ir 6 ee (4) yea (5.1) 


uw=f. y=ctgf(C,—C.f)+C,. 


Aa 2 ( an Peaks Or (5.2) 
oe fs y=tef (Cc, + C.f)+C,. 
y" + Fgxy'—2y=0 [Vegl. Kamxe (2.64)]. (5.3) 
Sonderfall von (5.1) mit w = arcsin (1/Coj x) . 
y = Sin x (C, — C, arcsin (1/Coj x)) + Cy. 
(a? — x2)?" = x (a? — x) y'—-2ary=0. (5.4) 
Sonderfall von (5.2) mit « = arcsin (%/a) . 
y = x[(a? — 37)!” (C, + C, arcsin(x/a)) + C,. 
a aye — se (1 = x4) — 2 Vim 0, (5.2) 


Sonderfall von (5.4) mit a=1. y= /(1— 2°)'? (C, + C, arcsinx) + Cy. 
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(1+ e-29) yy 2 y= 0. (5.6) 
Sonderfall von (5.2) mit « = arctg(e%* —1)”. 
y = (e2* —1)"? (C, + C, arctg (e?* — 11?) +. Cy. 
(x? +1) y"+2xy'—-2y=0. (5.7) 


Sonderfall von (5.2) mit w= arctgx, y= x (C,+ C,arctg x) +C,. 


Schlussbemerkungen 
Da v** F(u) — F**(u) v = 0 mit der Losung 


v= Fw) [C+ Cate du| 


die allgemeinste Form der reduzierten Differentialgleichung (3) darstellt, so ist 
die Differentialgleichung 


Fi) [»" - & y'] - 8 pry=0,. f=i@, f+0, ~ 
mit 
u=f und y=F() [C+ Cf (apy) | (6) 


die allgemeinste Differentialgleichung (1), die nach dem angefiihrten Verfahren 
losbar ist, sofern die beliebigen Funktionen F(f) und f(x) den eingangs gestellten 
Bedingungen geniigen. 

Die Beziehungen (5) und (6) sind einander eineindeutig zugeordnet. 


Summary 


By the transformation (4%) = v(u), u = f{ exp (—/ G(x) dx) dx the differen- 


tial equation D = y” + G(x) y’ + H(x) y = 0 turns to T = (u’)2 v** + Hv = 0, 


where v** signifies d?v/du®, and u’ = du/dw and H = H(x) should be expressed 
as functions of uw. 


From the solution v(“) of T follows immediately the solution y(*) of D, and 
vice versa. 


In this paper there are treated some of the types of differential equations, 
that may be solved by this method. 


(Eingegangen: 10. November 1954.) 
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On a Thermistor-Pirani Gauge of High Sensitivity 


By Ati Ati Arara, Fatut SuLTAN AHMED SULTAN, and IBRAHIM SAVED SHaFIE, 
Alexandria, Egypt?) 


In this paper an attempt is made to increase the sensitivity of a Pirani gauge 
by the introduction of a semi-conducting element having a negative temperature- 
resistance coefficient. The semi-conducting element which fulfills this requirement 
is a bead thermistor of the type D 163903 (V-519), used in 10-cm-band wave- 
guides for microwave-power detection. The electrical setup is shown in Figure 1, 


Figure 1 


in which 7, and 7, are 1000 and 250 2 respectively, 7, is the thermistor and 7, is a 
110 V, 25 W lamp which is connected to the vacuum system. It will be seen that 
the use of the thermistor in this connection has a great advantage as regards 
the gauge sensitivity and linearity. 

The direct-current characteristics of the thermistor are shown in Figure 2. 
As the current through the thermistor is increased, the voltage at first rises in 
accordance with Oum’s law. For large currents, however, appreciable heating 
of the thermistor results, thereby causing a decrease in resistance, and the cur- 
rent-voltage characteristic becomes non-linear. A sufficiently high current in- 
creases the temperature so much that the resulting drop in resistance may 
actually produce a voltage decrease and a negative differential resistance occurs. 

The best working conditions necessitate that the thermistor should be safe 
from overload or burnout. The useful and safe thermistor current ranges from 
about 2:5 to 11-0 mA, which corresponds to a resistance change from 670 to 1102 
and a voltage across it ranging from 1-58 to 1-2 V. 

When the pressure in the gauge increases, the resistance of the manometric 
wire decreases producing a drop in the potential of B relative to that of A. The 

current through the thermistor decreases producing an increase in its resistance. 
The potential of A accordingly becomes higher than it was, had the thermistor 
been replaced by a pure resistance. In other words, better sensitivity is obtained 
due to simultaneous changes of opposite signs in v, and 7, with changing pressure. 


1) Physics Dept., Faculty of Science, Alexandria University. 
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Figure 2 
Voltage-current and resistance-current characteristics of the bead thermistor. 


Figure 3 
(P—I)-relation for different voltages applied to the bridge; 


with thermistor and ——— with a pure resistance in arm 3. 
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Figure 4 


The relation between the pressure P and the thermistor current 7,, for different voltages applied 
to the bridge. 


500 


0 40 Pin 10mm Hg 80 
Figure 5 


The simultaneous variation of the filament and thermistor resistances with the pressure. 
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Figure 6 


The variation of the filament current with the pressure. 


The family of curves in Figure 3 demonstrate the relation between the pressure 
P in the gauge and the unbalanced current J in the meter M for different voltages 
(V) applied to the bridge. For each curve the thermistor current starts with a 
certain value, which is limited by the voltage V and the initial minimum pressure. 
These conditions fix the initial point on the working range BC of the thermistor- 
current-voltage characteristic. As the pressure increases, the thermistor current 
decreases (Figure 4). The working point moves towards the point B on the curve 
in Figure 2. The thermistor current values beyond B produced by further increase 
of pressure result in a non-linearity in the (P—J)-relation. V has to be so chosen 
that the initial thermistor current is high enough to realize a maximum linear 
range in the (P—J)-curve subject to the condition that this thermistor current 
is less than the burnout value. The relation is linear so long as the thermistor 
functions along the part BC. 
This corresponds to a thermistor resistance which ranges from 110 up to 670 Q. 
Figure 5 shows the variation of 7,, and rv, with P. The (P—v7,,)-relation is 
obtained by inserting a 1-Q-resistance milliameter in series with the thermistor 
and noting the variation of 7,, with increasing P, and then deducing the correspon-. 
ding values of 7, from the thermistor characteristics in Figure 2. The (P—7,)- 
relation, however, can be easily obtained by indirect means. 
Figure 6 gives the relation between P and i;. It is linear throughout that por- 
tion corresponding to the linear portion of the (P—I)-curve of Figure 3 and of th 
(P—1,,)-curve of Figure 4. 
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Sensitivity of the bridge 


In order to show the advantage realized by the thermistor, a similar set of 
curves (Figure 3) gives the relation between P and I using a pure resistance of the 
ame initial value as that of the thermistor at P,,,,,. They show a great differ- 
ence in the sensitivity and a wider range of linearity. The sensitivity of the gauge 
ising the thermistor amounts to 170 wA/ (10-* mm Hg), i.e. 5-88 X10-§ mm Hg/wA. 
rhe linearity of the curves surpasses that of the ordinary Pirani gauge. The sim- 
dlicity of the circuit compared with that of the ionization gauge is obvious. 

Such a bridge could be used with the utmost safety without the danger of burn- 
out for measuring the lowest pressure attained by any vacuum system, since the 
burnout current of the thermistor used is of the order of 50 mA. In the present 
work, a thermistor current of about 10-5 mA has not been exceeded. 

The choice of the pressure range is arbitrary, since by proper adjustment of 
R, a balanced circuit can always be obtained at any desired pressure, From 
this value of pressure, J increases linearly with P. This fact is indicated by the 
set of straight lines in Figure 7. When full scale deflection in the meter is reached, 
the balance is restored by decreasing the variable resistance FR in order to increase 
the current in the bridge and both filament current 7, and thermistor current 74, 
increase. 

The resistance of the filament decreases because the effect of increasing 
pressure overweighs that of increasing i;. The thermistor resistance decreases 


Pin 10 mm Hg (V=22V) 


40 
30 


20+ Cy 5) 


0 200 ; 400 microamperes 


Figure 7 
The (P—1)-relation for different values of R, 
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because of the increase in its current. Thus we start again with the same ratio 
between 7, and 7,. This process is repeated to cover the whole range of pressure 
required. With increasing P, the thermistor current decreases along one and the 
same (P—JI)-line. As mentioned above this has the effect of an augmented rise 
of the potential of A with respect to that of B. This set of lines are to a great 
extent parallel to each other, which corresponds to a constant sensitivity 
(1 wA/5-8 X10-§ mm Hg). Working from a pressure of the order of 2-7 X10—-? mm Hg 
downwards is possible using the same meter M (500 wA, 75 2 resistance), assum- 
ing the same sensitivity to apply to the extrapolated line OX. This assumption 
is justified, since the thermistor current is never below 3-5 mA, which is in the 
useful working range of the thermistor (Figures 2, 3, and 4). The minimu 
reading of the meter M in this case will depend on the efficiency of the vacuu 
system. 


Zusammenfassung 


Ein Hitzdrahtmanometer bildet, mit einem Thermistor in eine Wheatstone- 
Anordnung geschaltet, ein Thermistor-Pirani-Messgerat hoher Empfindlichkeit, 
das in einem weiten Gebiet eine lineare Druck-Strom-Abhangigkeit zeigt. Dabei 
wird von der verstarkenden Wirkung der negativen Widerstandcharakteristi 
eines Thermistors Gebrauch gemacht, um die von winzigen Druckaénderungen 
herriihrenden kleinen Widerstandvariationen zu verstarken. 


(Received: September 24, 1955.) 


The Indentation of a Transversely Isotropic Half-Space 
by a Rigid Punch 


By Harry Donatp Conway, Ithaca, New York, U.S.A.) 


Introduction 


Problems involving a large elastic mass having a plane boundary which i 
indented by a cylindrical punch whose axis is normal to the plane, have con- 
siderable practical significance. For example, they are encountered in the 
analysis of structural foundations where the large elastic mass simulates the 
earth, and the cylinder simulates a pillar. 

The first investigation of problems of this kind was made in 1885 by Bous- 
SINESQ [1]?) who gave the contact pressure distribution for a rigid, flat-ended and 


» Cornell University, College of Engineering, Department of Mechanics. 
) Numbers in brackets refer to References, page 84. 
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umber of punch shapes. In addition, ScHUBERT gave a very simple solution to 
he general problem of an indented isotropic half-plane and again illustrated his 
olution by a number of examples. 

Apparently unaware of SCHUBERT’s work, HarpING and SNEDDON [4] also 
reated the axisymmetrical problem, reduced it to that of solving dual integral 
quations different from those of ScHUBERT, and gave, as examples, the solutions 
or a flat-ended punch, cone, and ball. SNeppon later presented detailed numerical 
esults obtained by his method, including a study of the stresses in the interior 
if the half-space for a flat-ended punch [5] and cone [6]. An alternative method 
f solution, not confined to axisymmetrical problems, has been given by GREEN [7] 
vho has discussed, by way of illustration, the rigid cylindrical punch having a 
lat end inclined to the generators of the cylinder. 

Of recent years, a considerable amount of attention has been given to the 
yroblem of an indented half-space of transversely isotropic material, a number 
i papers having been published on the subject by Extiort [8], [9], SHiELp [10], 
ind PAYNE [11]. Such a material has the structure of hexagonal closely-packed 
rystals and is characterized by having five elastic constants. It is the purpose 
yf the present paper to show that, if the distribution of the contact pressure is 
cnown under an axisymmetrical rigid punch acting on an isotropic half-space, 
the corresponding distribution for a transversely isotropic half-space may be 
written down by inspection. Similar results [12] have recently been given for 
in indented orthotropic half-plane. 

As in most analyses of this kind, the effects of frictional forces under the 
junch have been neglected, and the plane assumed to be deformed solely by 
aormally applied forces. A two-dimensional investigation by OkuBo [13] indi- 
sates that this is quite probably a reasonable assumption for the axisymmetrical 
sase. It is also considered unlikely that the solutions are at all valid for an indent- 
ing surface of the punch which is far removed from a plane. For example, it is 
lifficult to see how it can reasonably be assumed that the applied forces are 
normal in the case of a cone with a small apex angle. Additionally, the problem is 
ussumed to be an elastic one, it being assumed that the presence of a small 
umount of plasticity will not change the calculated contact-pressure distribution 
to any great extent. 


Analysis 


Consider an isotropic half-space, the bounding plane of which is loaded with 
4 normal concentrated load P, acting in the z-direction. The solution of this 
problem is very well known, having first been given by BoussiNEsgQ [1], who 
showed that the vertical deflection of the bounding plane at a radius r¢ is 


Sao colt 2 a. af (1) 
neE v RE YY 


Og = 


with the usual notation for elastic constants. 

The case of a transversely isotropic half-space loaded in the same manner and 
having the bounding plane parallel to the hexagonal planes, was later treated 
by MicuELt [14], the corresponding vertical deflection being 


-{VACH LP (F+ 49" P 


2 
2nVL ACs Fs v (2) 


ZAMP VII/6 
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using the notation of Love [15] for the elastic constants A,C,F, and L, the strain 
energy function being given, in terms of these constants and the strains, by 
2W=A (62, + €2,) + C Chg te 2 \eyy ae) 
+ 2(A = 2.N) Cpptijg ch (02h Ce eae 
It follows that the vertical ‘deflection of the bounding plane in the transversely 


isotropic case may be immediately obtained from that in the isotropic case by 
merely replacing E’ by an equivalent modulus E,, defined by the equation 


(3) 


E 1 2 VA{WVAC as L)? — (P= yay a 
a2 ta 2 VL (AC ae F?) . 


Since the deformation caused by a rigid punch may be obtained by integration 
from the concentrated load solution if the contact pressure distribution is known 
it also follows that the pressure distribution for the transversely isotropic case is 
obtainable from that for the isotropic case merely by replacing E’ by E,,. 

In what follows, it will be advantageous to outline the solution of ScHu- 
BERT.[3], remembering now that the half-space may be tranversely isotropic 
rather than isotropic. Figure 1 shows the circular area of indentation caused by 


Figure 1 
Indented area of half space. 


an axisymmetrical punch. Denoting the pressure by p, the deflection at a radius 7 
due to the loading on the segments 4 BC and ADE is, from equation (2) 


, 


Ry 
Al 
to), == ae 1 
Wo nE,, dp [oR y) dR = x E,, dp ®,(¢) x (5) 
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Hence the deflection w, at a radius y due to the-€ntire loading on the circle is 


ee eat 
Oy RE a D,(p) dp. (6) 
0 
‘Now z = y sing and hence 
b.: 
| 2 D(z) mi ; ia OS 
a= SE / ae dz, @(z)= ®,(aresin =). (7) 


@(z) = 2 [ Pied ede | (8) 


It is interesting to note that the integral equation (8) has also been encountered 
by GuRNEE [16] in connection with optical studies on orientated filaments. The 
solutions of equations (7) and (8) are, respectively 


@O(z) =E es vy “_ @(") gg 
(2) = Ey, Sian rl, (9) 
a @ 
1 D(a) f @'(s) 1 Ga gO (Sins 
= os! FS Ry ey ee Ee 
P(e) a | Va? — o? J V/s? — @ ° mo de J Visti 0? : i 


2 


The constant term £,, w,(0) in equation (9) presents no difficulty, and is deter- 
mined from the equilibrium equation 


P=[22 (0) ode. (11) 
0 


Knowing the slope w{(v) of the punch, equation (9) gives ®(z) and hence Qs) 
which, when substituted in equation (10), gives the pressure distribution under 
the punch. 

Considering the flat-ended punch, w (rv) will be zero, ®(z) = @(s) = a constant, 
and the well-known result of BoussInEsQ [1] follows 


P(e) = x : (12) 


2a Va? — 9? 


It is seen that the indented area is constant, and the pressure is independent of 
the elastic constants of the half-space, being the same whether the latter is 
isotropic or transversely isotropic. This statement is only true in the case of a flat- 
ended punch, although the form of the contact pressure distribution will be the 
same for other kinds of punch. For example, the conical punch has a pressure 
distribution, determined from equations (9), (10), and (11) of 


Dk 


Ps arise prt all 13 
A eee cosh mab a aie cuca (13) 
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where « = tan, f being the angle which the face of the cone makes with the 
horizontal. 

It may also be mentioned here that BOROWICKA [17] has obtained, by a rather 
intuitive approach, the distribution of pressure under a flat-ended punch, when 
the load P is applied parallel to but distant b from the axis of the punch, Assum- 
ing that contact does not cease at any point, the contact pressure was shown 
to be % A 

0 
A(60) = = (1 +32. 8 cos 6] (14) 


2 wd Va? — 0 


where @ is measured from a diameter perpendicular to the axis about which the 
punch turns due to the eccentricity of loading. This result is also applicable to 
a transversely isotropic half-space. So far as is known, the solution has not been 
found for the case where the eccentricity b is sufficiently large to cause contact 
to cease over part of the face of the punch. However, it is apparent that the 
contact-pressure distribution is still the same whether the half-space is isotropic 
or transversely isotropic. 


REFERENCES 


[1] J. BousstnesQ, Application des potentials... (Paris 1885). 

[2] A.E.H.Love, Boussinesq’s Problem for a Rigid Cone, Quart. appl. Math. /0, 
161 (1939). 

[3] G.ScHuUBERT, Zur Frage dey Druckverteilung unter elastisch gelagerten Trag- 

__ werken, Ingen.-Arch. 73, 132 (1942). 

[4] J.W.Harpine and I.N.SNEppon, The Elastic Stvésses Produced by the 

Indentation of the Plane Surface of a Semi-Infinite Elastic Solid by a Rigid 

Punch, Proc. Camb. phil. Soc. 47, 16 (1945). 

[5] I.N.SnEppoN, Boussinesq’s Problem for a Flat-Ended Cylinder, Proc. Camb. 

phil. Soc. 42, 29 (1946). 

[6] I.N.SNEDDON, Boussinesq’s Problem for a Rigid Cone, Proc. Camb. phil. 

_ Soc. 44, 492 (1948). 

[7] A.E.GREEN, On Boussinesq’s Problem and Penny-Shaped Cracks, Proc. 
Camb. phil. Soc. 45, 251 (1949). . 

H. A. Evviort, Three-Dimensional Stress Distributions in Hexagonal Aeolo- 

0 tropic Crystals, Proc. Camb. phil. Soc. 44, 522 (1948). 

9} H.A.Ettiortt, A vial Symmetric Stress Distributions in Aeolotropic Hexagonal 

Crystals. The Problem of the Plane and Related Problems, Proc. Camb. phil. 

Soc. 45, 621 (1949). 

(20) Rea: SHIELD, Notes on Problems in Hexagonal Aeolotropic Materials, Proc. 

_ Camb. phil. Soc. 47, 401 (1951), 

[11] L.E. Payne, On Axially Symmetric Crack and Punch Problems for a Medium 
. with Transverse Isotropy, Proc. Camb. phil. Soc. 50, 466 (1954). : 

[12] Saas The Indentation of an Orthotropic Half-Plane, ZAMP 6, 402 

[13] H. OxuBo, The Stress Distribution in a Semi-Infinite Domain Having a 

Renae and Compressed by a Rigid Body, Z. angew. Math. Mech. 20, 

{14} J.H.Micnerr, The Stress in an Aeolotropic Elastic Solid with an Infinite 

Plane Boundary, Proc. Lond. math. Soc. 32, 247 (1900). 


[15] A. E.H. Love, The Mathematical Theory of Elasticity, 4th edition (Cambridge 
University Press 1927), 160 p. 


- 


Vol. VII, 1956 Varia — Miscellaneous — Divers 85 


[16] E.F.GuRNEE, Solution of an Integral Equation A rising in Optical Studies of 
Oriented Filaments, J. appl. Phys. 26, 918 (1955); 

[17] H. Borowicxa, Uber ausmittig belastete starve Platten auf elastisch-isotropem 
Untergrund, Ingen.-Arch. 74, 1 (1943) 


Zusammenfassung 


Der Verfasser gibt eine allgemeine Lésung fiir die Verteilung des Druckes 
zwischen einem axialsymmetrischen Stempel und einem transversal-isotropen 
Halbraum. Es wird gezeigt, dass die Verteilung des Druckes fiir den flachen 
Stempel mit allgemeiner Belastung unabhangig ist von den elastischen Eigen- 
schaften des Halbraums und auch genau dieselbe, als ob der Halbraum isotrop 
ware. 


(Received: September 13, 1955.) 
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Congrés international sur les tubes électroniques pour hyperfréquences 
du 29 mai au 2 juin 1956 4 Paris 


La Société des radioélectriciens et la Société frangaise des Ingénieurs tech- 
niciens du vide organisent a Paris, du 29 mai au 2 juin 1956, le congrés sub- 
mentionné. Les sujets proposés sont: 


. Phénoménes fondamentaux et théories générales. 

. Lampes a grilles de contrdéle. 

. Magnétrons. 

Tubes 4 modulation de vitesse. 

Tubes 4 onde progressive amplificateurs et oscillateurs. 

. Tubes divers (tubes a saut de potentiel, tubes 4 ondes de charge d’espace...). 
. Accélérateurs linéaires. 

. Générateurs spéciaux pour ondes millimétriques et inframillimétriques. 
. Tubes a décharge gazeuse. 

10. Circuits associés (cavités résonnantes, circuits a retard). 

11. Technologie. 

12. Optique électronique. 

13. Mesures. 


Siége du Congrés: Paris VI°, 44, rue de Rennes. 
N. SCHAETTI 


9. Internationaler Kongress fiir Angewandte Mechanik 


Der Kongress wird vom 7. bis 13. September 1956 in Briissel stattfinden. 
Kongressleitung: Prof. Dr. F. vaN DEN DUNGEN, Universitat Briissel. Sektionen: 
1. Hydro- und Aerodynamik; 2. Mechanik der festen K6rper (starre Korper, 
Schwingungen, Elastizitat, Plastizitat). Das Programm kann von den Lehr- 
stiihlen fiir Mechanik an der ETH. bezogen werden. H. ZIEGLER 
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Primzahlen. Par Ernst Trost (Birkhauser, Bale 1953). Volume 2 de la série: 
Elemente der Mathematik von héherem Standpunkt aus. 95 p.; broché fr.s. 13.50. 

On sait que la théorie des nombres premiers pose des problémes intéressants, 
dont quelques-uns sont encore a résoudre. EUCLIDE a déja établi l’existence d’une 
infinité de nombres premiers; mais ce n’est qu’au XIX* siécle que DIRICHLET a 
démontré qu’une progression arithmétique, dont les termes sont premiers avec la 
raison, contient aussi une infinité de nombres premiers. A la fin du méme siécle 
HapAMarp et DE LA VALLEE-Poussin établirent la loi asymptotique de la distri- 
bution des nombres premiers: si 2(¥) désigne le nombre des nombres premiers 
inférieurs & x, le quotient de 2(*) par #/In x tend vers 1 lorsque # croft indé-— 
finiment. Une proposition célébre de GoLpBACcH ~ il existe une infinité de nombres 
premiers p tels que p+ 2 soit encore un nombre premier — n’est pas encore 
démontrée. 

L’opuscule de M. Trost est une excellente introduction a la théorie des nom- 
bres premiers et aux questions rappelées ci-dessus. En particulier la loi asymp- 
totique de la distribution des nombres premiers et le théoréme cité de DIRICHLET 
y sont établis par la méthode récente de SELBERG et d’ERD6s, qui constitue une 
simplification remarquable et est de nature élémentaire. L’exposition est concise 
et rigoureuse. Bien que demandant du lecteur une attention et un effort soutenus, 
la lecture ne réclame de lui que peu de connaissances préalables d’analyse. Les pro- 
priétés de théorie des nombres nécessaires sont établies dans la. mesure ot elles” 
trouvent emploi pour celle des nombres premiers. L’ouvrage peut étre recom- 
mandé aux jeunes mathématiciens qui désirent étre mis au courant des questions ~ 
les plus importantes de la théorie des nombres premiers. M. Plancherel 


Les fonctions orthogonales dans les problémes aux limites de la phy- 
sique mathématique. Von TH. VoGeEL (Editions du Centre national de la re- 
cherche scientifique, Paris 1953). 191 S., 4 Fig., broschiert fFr. 1200.—. 

Die ausserordentliche Bedeutung der orthogonalen Funktionssysteme fiir die 
Mathematik und ihre Anwendungen rechtfertigen wohl die Herausgabe eines 
Leitfadens tiber diesen Gegenstand. Die vorliegende Schrift behandelt zuerst 
die allgemeine Theorie unter Einschluss der Eigenwertprobleme, die ja auch auf 
orthogonale Funktionssysteme fiihren. In einem 2. Kapitel werden die klassi- 
schen Orthogonalsysteme (Besselsche Funktionen, Legendresche, Hermitesche, 
Laguerresche Polynome und andere) und insbesondere auch die Fourierschen 
Reihen eingehend behandelt. Das 3. Kapitel enthalt sodann die zahlreichen, sorg- 
faltig ausgewdhlten und durchaus nicht trivialen Anwendungsbeispiele aus der 
Physik. Das Werk kann zum Einarbeiten in die Materie warmstens empfohlen 


werden. H, Rutishauser 


Discontinuous Automatic Control. Von InmGarD FLiGGr-Lorz (Princeton 
University Press, Princeton, N. J., 1953). 168 S., 102 Fig.; $5.-. 

In vielen Fallen liesse sich eine Regulierung technisch am einfachsten lésen, 
wenn das Regelorgan nur zwei symmetrische Stellungen oder zwei entgegen- 


gesetzt gleiche Verstellgeschwindigkeiten aufweisen miisste. Eine solche Lésung © 
lasst sich in geschlossener mathematischer Form nicht mehr diskutieren. Im vor- — 


ett tei 
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liegenden Buch wird versucht, den Ausschwingvorgang, die Stabilitatsbedin- 
gungen und den Einfluss von Unvollkommenheiten des Regelsystems systema- 
tisch mit Hilfe graphischer Verfahren eingehend zu diskutieren. Drei Viertel des 
Buches sind der Betrachtung eines bewegten K6rpers mit einem Freiheitsgrad 
gewidmet, wobei die Aufgabe der Regelung darin besteht, beim Auftreten einer 
dussern Stérung den Kérper méglichst rasch in die urspriingliche Lage oder 
Bahn zuriickzufiihren. Die Stellung (Positionsregelung) bzw. die Verstellge- 
schwindigkeit (Geschwindigkeitsregelung) des Regelorgans ist nur vom Vor- 
zeichen einer Regelfunktion abhangig, welche ihrerseits von der Lageabweichung 
und deren zeitlicher Ableitung des betrachteten K6rpers abhangt. Die der Haupt- 
bewegung iiberlagerte Bewegung soll sich durch eine lineare Differentialgleichung 
zweiten Grades darstellen lassen. In diesem Falle lasst sich die superponierte 
Bewegung fiir jede Ausgangslage in zeitliche Teilstiicke zerlegen, wobei die Zeit 
als Parameter erscheint und die Bewegung in einem Phasendiagramm (Lageab- 
weichung und deren zeitliche Ableitung als schiefwinklige Koordinaten) als 
Stiicke einer logarithmischen Spirale aufgezeichnet werden kann. Die Diskonti- 
nuitat hat zur Folge, dass nur zwei Spiralen mit Zentren symmetrisch zum 
Ursprung betrachtet werden miissen. Die Schaltpunkte liegen auf einer Geraden 
durch den Ursprung des Phasendiagramms (Positionsregelung). Bei der Ge- 
schwindigkeitsregelung mit Riickfiihrung im Kontrollsystem wird gezeigt, dass 
die Darstellung in zwei Ebenen eines Phasenraumes geniigt, um alle Falle disku- 
tieren zu kénnen. In dieser Darstellungsart lasst sich der Einfluss von Unvoll- 
kommenheiten des Regelsystems (zeitliche Umschaltverzégerung; Umschaltung 
erst bei einem bestimmten Wert der Regelfunktion; oder Regelorgan auf Null, 
sofern die Regelfunktion innerhalb bestimmter, entgegengesetzter Werte sich 
befindet) sehr gut beriicksichtigen. Im letzten Teil des Buches wird dieselbe Art 
Betrachtung angewandt auf das Problem der Langsstabilisierung eines Flug- 
k6rpers mit vorgeschriebener gerader Bahn. Die Regelung geschieht durch 
rasche Verdrehung des Einstellwinkels des Tragfliigels oder des Héhenleitwerks 
zwischen zwei bestimmten Werten in Abhangigkeit der Abweichung von einem 
vorgegebenen absoluten Winkel und deren zeitlicher Ableitung. Die Bewegung 
selbst hat drei Freiheitsgrade; ohne Regelung ergeben sich vier Lésungen der 
freien Bewegung. Bei Annahme von je zwei konjugiert komplexen Lésungen 
wird gezeigt, wie die Bewegung mit der Regelung unter verschiedenen Bedin- 
gungen erfolgt. Ein numerisches Beispiel wird durchgerechnet und in Phasen- 
diagrammen diskutiert. 

Das Wertvolle an diesem Buch ist die sehr eingehende Darstellung der Me- 
thode, mit welcher man den Problemen mit diskontinuierlicher Regelung bei- 
kommen kann. Literaturhinweise fehlen fast vollkommen, was sich damit er- 
klaren liesse, dass sich sozusagen alle Arbeiten iiber dieses Gebiet hinter dem 
Vorhang der militarischen Geheimsphare verborgen halten. Um so dankbarer 
darf die zivile Welt iiber diesen Beitrag sein, der eine Liicke der Servotechnik 
ausfiillen hilft. H, Weber 


Gasdynamik. Von K. Oswartitscu (Springer-Verlag, Wien 1952). 4565., 
300 Abb.; Ganzleinen sFr. 80.—. 

Die Gasdynamik, die Lehre von den Strémungen kompressibler Medien, hat 
in den letzten zwei Jahrzehnten einen sehr bedeutenden Aufschwung genommen. 
Vor allem ist es die Flugtechnik, die, nachdem durch neuartige Antriebe viel 
hdhere Geschwindigkeiten méglich geworden sind, sich sehr eingehend mit der 
Kompressibilitat befassen muss. 


1 
{ , 
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Das Buch erscheint zu einer Zeit, in der das Bediirfnis nach Ordnung und 
Verstindlichmachung der sehr zahlreichen theoretischen und experimentellen 
Untersuchungen der Kriegszeit besonders gross ist. Da der Verfasser es ausge~ 
zeichnet versteht, mit anschaulichen Begriffen und Methoden zu arbeiten, handelt 
es sich um ein Werk, das von grésstem Wert fiir jeden ist, der sich in dieses neue 
Gebiet einarbeiten will. 4 

Daneben zeichnet es sich durch eine aussergewohnliche Vollstandigkeit der 
Erfassung der zum Teil schwer zugdnglichen Literatur aus. Fiir jeden Kenner 
der Lehrbuchliteratur ist es eine angenehme Uberraschung, zu sehen, wie OSWA~ 
tTitscH keineswegs nur die bekannten Uberlegungen wiedergibt, sondern nicht 
selten recht originelle Gedankenwendungen vornimmt, die den Gegenstand in 
neuem Licht erscheinen lassen. ; 

Nach einleitenden Kapiteln mehr physikalischer Natur wird zunachst die 
Fadenstromung behandelt, wobei die auch hier niitzlichen Charakteristiken friih- 
zeitig eingefiihrt werden. Insbesondere ist der Abschnitt iiber nichtstationare 
Fadenstrémung sehr ausfiihrlich, indem zum Beispiel die kugel- und zylinder-_ 
symmetrischen Strémungen mit und ohne Verdichtungsstésse behandelt werden. 
Sodann folgen die allgemeinen Satze iiber zwei- und dreidimensionale Stromun- 
gen. — Unter den Integralsdtzen ist der Oswatitschsche tiber die Widerstands- 
berechnung aus dem Entropiefluss besonders bemerkenswert. : 

Sehr ausfiihrlich wird sodann die linearisierte Theorie der Unter- und Uber- 
schallstrémungen dargestellt. Es werden transsonische Ahnlichkeitsgesetze be-— 
handelt, ebenso die Hyperschallstr6mungen (sehr grosse Mach-Zahlen). 

Ein grésserer Abschnitt ttber kompressible Grenzschichten und ein Kapitel 
iiber Messtechnik beschliessen das Buch. — Druck und Ausstattung sind in gewohn- 
ter Weise vorziiglich. Man darf dem Buch eine weite Verbreitung wiinschen. 

J. Ackeret 


Some Basic Problems of the Mathematical Theory of Elasticity. Von 
N. I. MUSKHELISHVILI; tibersetzt aus dem Russischen durch J. M. R. RaDoK 
(P. Noordhoff Ltd., Groningen, 1953). 704 S., 65 Fig.; hfl. 38.-. : 

Der Titel dieses Buches ist allzu bescheiden. In Wirklichkeit enthalt das Werk 
_ neben einer vorbildlichen Einfiihrung in die lineare Elastizitatstheorie mit An- 
wendungen auf die exakte Balkenlehre eine umfassende Darstellung der mathe- 
mathischen Theorie der ebenen Zusténde (Spannungs- und Verzerrungszustand) 
in isotropen K6rpern auf funktionentheoretischer Basis. P- 

Die komplexe Behandlung ebener Zustande geht auf G. V. KoLosow zuriick 
und ist vorwiegend von russischen Forschern zu einer Theorie ausgebaut worden, 
die der komplexen Strémungslehre zur Seite zu stellen, mit Riicksicht auf das 
differenziertere Verhalten des elastischen Kérpers, aber wesentlich verwickelter 
ist. Das Buch enthalt neben einer neuen Darstellung der Kolosowschen Theorie 
eine bemerkenswerte Fiille von neuen Anwendungen allgemeiner wie spezieller 
Natur, die sich aus der systematischen Verwendung der Eigenschaften der 
Cauchy-Integrale sowie der Méglichkeiten der konformen Abbildung ergeben. 
Besonders sorgfaltig sind die Existenz- und Eindeutigkeitsfragen behandelt, die 
bei mehrfach zusammenhangenden Bereichen (Ringe) wichtig sind. ey 
Nees ise vorauszusehen, dass dieses Werk auf Jahre hinaus auf seinem Gebiete 
grundlegend bleiben wird. Mit seiner sorgfaltigen (ubrigens auch buchtechnisch _ 
vorbildlich geratenen) Ubertragung in die englische Sprache hat sich auch der 
Ubersetzer ein grosses Verdienst erworben. ' -H. Ziegler 
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Aus dem Inhalt: Uber gewohnliche Differentialgleichungen mit Singularititen und die zuge- 
h6rigen Entwicklungen willkiirlicher Funktionen — Das asymptotische Verteilungsgesetz der 
Bigenschwingungen eines beliebig gestalteten elastischen Kérpers — Uber die Gleichverteilung 
yon Zahlen mod. Eins — Uber die Bestimmung einer geschlossenen konvexen Fliche durch 
ihr Linienelement — Gravitation und Elektrizitit —- Ausbreitung elektromagnetischer Wellen tiber 
einen ebenen Leiter — Uber die neue Grundlagenkrise der Mathematik — Zur Infinitesimalgeo- 
mettie: Hinordnung der projektiven und der konformen Auffassung — Theorie der Darstellung 
kontinuierlicher halb-einfacher Gruppen durch lineare Transformationen (Teil I, II, III und 
Nachttag) — Integralgleichungen und fast-periodische Funktionen — Die Vollstaindigkeit der 
primitiven Darstellungen einer geschlossenen kontinuierlichen Gruppe (mit F. Peter) — Uber das 
Pick-Nevanlinnasche Interpolationsproblem und sein infinitesimales Analogon — Spinots in 
Dimensions (with R. Brauer) — The Method of Orthogonal Projection in Potential Theory — 
On the Differential Equations of the Simplest Boundary-Layer Problems — On Hodge’s Theory 
of Harmonic Integrals — Fundamental Domains for Lattice Groups in Division Algebra (Part I 
and II) — Elementary Algebraic Treatment of the Quantum Mechanical Symmetry Problem — 
Die natiitlichen Randwertaufgaben im Aussenraum fiir Strahlungsfelder beliebiger Dimension. 
und beliebigen Ranges — Verzeichnis der Publikationen Hermann Weyls bis 1955. 


Diese Auswahl der bedeutendsten Abhandlungen H. Weyls aus der Zeit von 1910 bis 1952 gibt 
einen Querschnitt durch das mathematische Werk des grossen Gelehrten und widerspiegelt die 
Entwicklung fast aller Teile der Mathematik wahrend einer langen und wichtigen Periode. 
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